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ABSTRACT:  Jklre-^laTlSTalmed-  to--b»-the  first  -con a&e&flT  attempt.To  match 
experimental  and  theoretical  data  on  long-distance  propagation  of 
superlong  (wave-lengths  of  several  times  ten  kilometers)  waves  around 
the  earth.  Since  there  is  no  probability  distribution  function  for 


the  experimental  data,  the  matching  is  carried  out  approximately  by  the 
method  of  mixed  initial  data,  where  all  the  data  on  the  field  and  on 
the  medium  are  divided  into  two  groups  -  reliable  and  unreliable. 

Several  models  are  proposed  for  the  propagation  along  the  earth's  surface 
and  in  the  ionosphere  with  an  attempt  to  include  all  the  geophysical 
factors  which  influence  the  far  field  of  superlong  radio  waves.  Only  the 
waveguide  channel  adjacent  to  the  earth  is  considered.  The  method  of 
normal  waves  which  can  be  used  to  solve  waveguide  propagation  problems 
for  sound  waves  in  the  ocean,  Infrasound  waves  in  the  atmosphere,  seismic 
waves  in  the  earth,  etc.  is  also  developed. fJ  The  first  edition  was 
published  in  1956.  English  translation:  ltfj.  pages. 
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PREFACE  TO  SECOND  EDITION 
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1.  More  than  seven  years  have  passed  since  first  publication  of  this  book, 

and  it  went  out  of  print  in  1956.  In  order  to  decide  question  about  expediency  of 

* 

its  republication  one  should  consider  material  contained  in  light  of  new  works  on 
propagation  of  very  long  radio  waves  and  related  geophysical  problems  and  should 
also  examine  changes  which  leading  physical  concepts  and  mathematical  methods  that 
are  basis  of  this  book  have  undergone  during  the  last  few  years. 

2.  From  contemporary  point  of  view  [1]  problem  of  theoretician  studying 
propagation  of  radio  waves  under  terrestrial  conditions  consists  in  coordination  of 
experimental  data  on  field  of  waves  and  geophysical  data  essentially  affecting  radio 
waves  with  help  of  functional  relationships  between  these  data.  In  order  that 
theory  be  useful  for  practical  workers,  functional  relationships  between  field  and 
medium  have  to.be  derived  from  equations  of  field  and  material  equations  of  medium 
and  not  directly  from  experimental  data.  These  relationships  establish  one-to-one 
correspondence  between  functions  characterizing  medium  and  functions  describing 
electromagnetic  fields.  In  ideal  case,  when  experimental  data  are  in  the  form  of 
distribution  function  of  probability,  determining  degree  of  authenticity  of 
experimental  data,  coordination  consists  in  imposition  on  arguments  of  this  function 
of  theoretical  functional  relationship,  which  transforms  unconditional  probability 
into  conditional.  Thus  is  obtained  gain  of  information  according  to  Shannon,  since 
initial  regions  of  uncertainty  of  experimental  data  narrow.  Similar  state  of  affairs 
exists  in  statistical  theory  of  information  transmission  when  a  priori  knowledge  of 
information  is  definltlzed  according  to  data  on  signal  carrying  the  Information  [2]. 
Let  us  note  that  in  process  of  coordination  of  data  there  appears,  gain  of  information 
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both  for  field  and  for  medium.  Its  magnitude  depends  on  character  of  functional 
relationship  and  degree  of  authenticity  of  experimental  data.  As  it  is  known, 
information  may  be  obtained  only  from  experiment.  In  this  case  no  additional 
experiments  are  carried  out  and  gain  of  information  is  due  to  experimental  data 
Included  in  theoretical  functional  relationships.  Progress  of  theory  consists  in 
that  owing  to  discovery  of  new  facts,  more  precise  definition  of  old,  or  desire  to 
include  known  facts  in  realm  of  coordination,  there  occurs  from  time  to  time 
recoordination  of  data  on  field  and  medium  on  the  basis  of  more  exact  or  broader 
functional  relationships. 

This  book,  as  far  as  I  know,  constitutes  first  conscious  attempt  to  apply 
concepts  of  coordination  to  problem  of  long-range  radio  communications  on  kilometer 
waves  (very  long)  around  Earth.  Owing  to  absence  of  distribution  functions  of 
probability  for  experimental  data,  coordination  was  produced  by  approximation  method 
of  mixed  initial  data.  It  consists  in  separation  of  all  data  on  field  and  medium 
into  two  groups:  reliable  and  unauthentlc  data.  First  are  introduced  into 
functional  relationships  as  initial  data,  while  second  are  determined  from  these 
relationships  as  unknown  quantities.  Method  of  mixed  initial  data  gives  certain 
loss  of  information  as  compared  to  ideal  method  of  coordination,  since  in  it 
information  contained  in  unlikely  data  is  disregarded  and  reliability  of  reliable 
data  introduced  as  initial  is  reevaluat'd.  However,  this  loss  is  not  so  big  as 
obtained  during  solution  of  primal  problems  of  propagation  of  redio  waves  for 
which  data  on  field  are  generally  disregarded  and  all  data  on  medium  are  considered 

f 

absolutely  reliable.  For  problem  of  propagation  of  very  long  radio  waves  around 
Earth  there  is  obtained  in  this  case  negative  value  for  gain  of  information  according 
to  Shannon  [1]. 

5.  Data  obtained  as  a  result  of  coordination  will  have  meaning  of  new 
information  on  field  and  medium  only  under  the  condition  that  functional  relationships 
are  adequate  to  true  ones.  There  was  possibility  of  obtaining  uncontradictory 
coordination  of  limited  number  of  experimental  data  says  nothing  for  adequacy  of 
functional  relationships.  Coordination  can  also  be  obtained  with  Incorrect 
functional  relationship,  if  it  contains  sufficient  number  of  arbitrary  constants. 

As  exavole  of  this  we  have  numerous  phenomenological  formulas  for  calculation  of 
distant  fi?ld  of  radio  waves  (Osteen,  Ekspenshid,  and  others)  appearing  till  now  in 
many  textbooks  and  monographs  on  propagation  of  radio  waves.  Another  method,  creating 
illusion  of  adequacy  of  functional  relationships,  consists  in  that  some  of  parameters 
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of  medium  or  field  are  declared  effective,  and  it  is  permitted  to  take  different 
values,  depending  upon  what  group  of  facts  is  subjected  to  coordination.  Such 
method  was  used  during  the  study  of  lower  ionosphere  [9-H]  (points  4  and  5  this 
preface),  which  in  recent  past  was  still  in  the  full  sense  a  blank  spot  in 
geophysics,  and  it  was  possible  to  ascribe  any  properties  to  it.  Necessity  to 
change  parameters  of  ionosphere  during  transition  fron  one  particular  coordination 
of  data  to  another  is  criterion  of  inadequacy  of  model  of  medium.  This  criterion 
was  used  in  present  work  during  construction  of  models  A  and  B  (Table  1) .  However, 
basic  method  of  selection  of  correct  model  of  medium,  requiring  enlistment  of 
minimum  quantity  of  experimental  data,  consisted  in  the  following.  Having  assigned 
a  group  of  facts  to  be  subjected  to  coordination  and  having  selected  certain  model 
of  medium,  we  estimated  what  parameters  of  model  essentially  influence  these  facts 
through  functional  relationships.  Gradually  complicating  model  of  medium  and 
constructing  functional  relationships  between  field  and  medium  for  it,  we  developed 
simplest  models  of  medium  A  and  B  (Table  1),  considering  all  basic  geophysical 
factors  essentially  affecting  properties  of  distant  field  of  very  long  radio  waves. 

As  a  result  of  such  analysis  we  determined  necessity  for  calculation  of  sphericity 
of  earth  and  ionosphere,  conductivity  of  earth,  magnetic  anisotropy  of  ionosphere, 
and  also  changeability  of  medium  along  route  of  radio  wave.  Lower  part  of  ionosphere 
was  replaced  by  step  model,  characterizing  concentration  of  electrons  and  effective 
frequency  of  collisions  of  electrons  with  neutral  molecules.  Possibility  of 
disregarding  motion  of  ions  in  lower  ionosphere,  in  spite  of  their  large  quantity 
in  this  region,  follows  from  works  of  Bates,  Messi,  and  Nicolet.  Effective 
frequency  of  collisions  was  determined  from  Nicolet' s  theory,  which,  as  recent 
investigations  showed  [5],  gives  somewhat  exaggerated  values  in  connection  with 
fact  that  transport  cross  section  for  molecules  of  nitrogen  is  proportional  to 
velocity  of  electrons  and  not  inversely  proportional  to  it  as  was  assumed  by  Nicolet. 
Doubt  can  also  be  caused  by  replacement  of  real  smooth  distributions  of  electron 
concentration  and  frequency  of  collisions  with  height  by  one-  or  two-stage  curves. 
However,  more  general  theory  of  propagation  of  long,  very  long,  and  extremely  long 
radio  waves  around  Earth  [1,  4]  showed  that  such  idealization  xs  permissible  for 
waves  with  frequency  f  <  JO  kilocycles  and  distances  D  >  500-1000  km.  In  these 
works  it  was  shown  that  for  expansion  of  functional  relationship  to  region  of  data 
embracing  frequency  up  to  100  kc  and  any  distances,  introduction  of  several  steps 
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is  inexpedient.  More  effective  results  are  given  by  smooth  model  of  electron 
concentration  and  frequency  of  collisions.  In  these  works  was  also  clarified 
unreasonableness  of  use  as  initial  data  of  measurements  of  near  field  only.  This 
is  connected  with  fact  that  on  long-range  routes  there  can  be  considerable  local 
heterogeneities,  caused  by  geoactivity  of  Sun  and  radiation  belts  of  Earth.  Thus 
appraisal  of  possibility  of  predictions  of  properties  of  distant  field  according  to 
data  on  near  field  only  is  too  optimistic. 

4.  Consideration  of  problem  in  light  of  information  theory  has  obscured  the 
semantic  side  of  question.  At  the  same  time,  as  a  result  of  coordination  with  help 
of  adequate  functional  relationships  between  field  and  medium,  the  mechanism  of 
long-range  propagation  of  very  long  radio  waves  around  Earth  is  clarified.  As 
follows  from  Chapter  ?,  distant  field  is  in  the  form  of  spectrum  of  normal  waves 
[5j  8]  propagating  in  unique  waveguide  shaped  like  surface  of  earth  and  lower  layers 
of  ionosphere.  These  waves  are  similar  to  TH  and  TE  waves  in  hollow  metallic  pipes, 
known  since  times  of  Rayleigh  [6],  However,  Watson,  having  obtaining  in  1919  [7]  a 
field  of  long  waves  in  the  form  of  spectrum  of  normal  waves,  did  not  discover 
relationship  of  propagation  in  hollow  pipes  to  waveguide  mechanism.  Method  of 
normal  waves  [5,  8]  utilized  in  this  book  establishes  this  relationship  automatically. 
Therefore  we  use  waveguide  terminology  throughout.  Thanks  to  existence  of  critical 
frequencies  for  all  waves  but  THq,  determined  by  ratio  of  width  of  waveguide  to 
wavelength,  over  large  distance  of  field  of  radio  wave,  there  are  not  more  than 
two-three  normal  waves.  It  is  necessary  to  note  that  under  real  conditions,  besides 
surface  waveguide,  there  are  underground  waveguide,  through  which  travel  normal 
waves,  creating  above  earth  so-called  lateral  waves  (terminology  of  Ott  and  L.  M. 
Breknovskikh) ,  and  ionospheric  waveguides,  among  which  is  separated  waveguide 
channel  of  whistling  atmospherics  [sic].  General  diagram  of  spectra  of  normal  waves 
is  analyzed  in  works  [1,  4],  In  this  book  only  surface  waveguide  channel  is 
considered. 

We  direct  attention  of  reader  to  necessity  for  calculation  of  magnetic  field 
of  Earth,  sphericity  of  waveguide  walls,  and  finite  conductivity  lower  wall  of 
Earth.  Magnetic  anisotropy  of  ionosphere,  appearing  in  night,  cancels  Brewster 
effect.  Without  it  distant  field  of  very  long  waves  would  be  weaker  at  night.. 
Disregard  for  sphericity  of  Earth  and  ionosphere  at  very  long  lengths,  and  all  the 
more  so  at  long  wavelengths,  is  absolutely  impermissible.  This  is  connected  with 
existence  of  effect  of  adhesion  of  waves  of  any  nature  to  concave  surfaces:  in  this 
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case,  to  lower  ionosphere.  Effect  of  adhesion  was  discovered  with  acoustic  waves 
by  Rayleigh;  with  radio  waves  it  was  investigated  in  work  [5,  8],  Normal  waves 
in  which  effect  of  adhesion  appears  have  angular  wave  numbers  lying  in  interval 
between  ka  and  kc,  c  «  a  +  h;  a  -  63 70  km.  This  means  that  they  possess  phase 
velocities  less  than  velocity  of  light,  as  can  be  seen  from  Pig.  17.  In  waveguide 
with  flat  walls  phase  velocities  of  normal  waves  are  always  greater  than  velocity  of 

light,  unless  we  ascribe  supernatural  properties  to  earth  or  ionosphere.  Adhesion 

* 

of  normal  waves  to  ionosphere  considerably  weakens  coefficients  of  excitation  of 
normal  waves  nj  (see  Fig.  19)  as  compared  to  flat  case.  Comparisons  of  flat  and 
spherical  waveguide  are  given  in  Pig.  3  of  work  [1].  With  increase  of  frequency, 
besides  effect  of  adhesion,  there  occurs  weakening  of  reflectance  of  ionosphere, 
and,  as  a  result  of  these  two  effects,  mechanism  of  propagation  of  waves  in 
continuous  manner  passes  from  waveguide  to  diffraction.  This  transition  is 
accompanied  by  replacement  of  normal  waves,  dominating  at  large  distance  from 
radiator.  Thus,  for  instance,  for  extremely  long  waves  (f  <  3000  cps)  in  distant 
field  THq  wave  dominates,  also  called  TEM  or  cabin  wave.  In  range  of  frequencies 
from  10  to  20  kilocycles,  considered  in  this  book,  there  dominates,  as  a  rule,  the 
TT^  wave,  while  in  certain  cases  the  TH2  is  dominant.  At  frequencies  of  50-60 
kilocycles  by  day  TH2  and  TH^  waves  dominate,  and  at  night  the  TH^  and  TH^  are 
dominant  [1,  4].  As  a  result  of  replacement  of  leading  waves  average  phase  velocity 
along  Earth  at  large  distances,  determined  by  leading  wave,  is  always  greater  than 
velocity  of  light  for  f  >  1000-2000  cps.  Diurnal  variations  of  distant  field  of 
vecy  long  and  long  waves  are  determined  by  interference  and  replacement  of  above- 
indicated  normal  waves.  Finally,  in  present  work  is  established  influence  of  finite 
conductivity  of  earth  on  attenuation  factors  of  leading  normal  waves  (see  Fig.  15). 
For  "soil"  of  average  conductivity  they  are  increased  by  approximately  0.5  neper  per 
radian  of  distance. 

5.  However,  attempts  continue  to  insert  certain  experimental  data  on  very 
long  waves  into  the  Procrustean  bed  of  functional  relationships  derived  from 
simplified  models  of  medium  in  which  are  disregarded  sphericity  of  earth  and 
ionosphere,  conductivity  of  earth,  magnetic  anisotropy  of  ionosphere,  and 
changeability  of  medium  along  route  of  wave.  Impossibility  of  uncontradictory 
coordination  with  help  of  these  functional  relationships  follows  in  certain  cases 
from  works  themselves.  Thus,  for  instance,  in  works  (9,  10,  11]  the  author  had  to 
change  height  of  lower  layer  of  ionosphere  from  70  to  50  km  and  its  conductivity  in 
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range  10  -  10  CGSF,  depending  upon  selected  frequency  and  distance  of  route,  which, 

as  noted  above,  Is  first  criterion  of  Inadequacy  of  model  and  was  used  by  us  for  its 
correction.  Thus  formulas  obtained  by  them  differ  little  from  phenomenological 
formulas,  thanks  to  contribution  of  Osteen,  Ekspenshid,  and  others.  The  route  to 
such  simplification  was  taken  also  by  certain  foreign  authors  [20  and  21],  Only 
recently  have  works  considering  some  of  esstential  factors  separately  began  to 
appear  [12,  13],  From  Chapters  IV  and  V  of  this  book  it  follows  that  these  i actors 
interact  with  each  other  and  that  their  calculation  separately  cannot  bring  benefit 
to  practice.  Thus,  for  instance,  calculation  of  conductivity  of  earth  in  flat  model 
yields  exaggerated  attenuation  factors  [20],  since  effect  of  adhesion  weakens 
influence  of  earth. 

6.  From  the  point  of  view  of  contemporary  state  of  problem  of  propagation  of 
very  long  radio  waves  around  Earth,  in  this  book  the  following  results  are  the  most 
important:  1)  explanation  of  diurnal  variations  of  distant  field  of  very  long 
waves;  2)  prediction  of  diurnal  variations  in  phase  of  distant  field,  confirmed 
completely  by  experiments  of  Pierce  and  Heritage  in  the  United  States  [14,  15],  and 
selection  of  the  most  phase-stable  range  of  waves  useful  for  transmission  of  signals 
of  exact  time;  3)  appraisal  of  electron  concentration  of  lower  ionosphere  by  day  and 
at  night.  With  respect  to  the  latter,  oue  should  note  that  replacement  of  smooth 
distribution  of  profile  of  electron  concentration  by  stepped  type  permits  us  to 
treat  equivalent  parameters  only.  In  reality,  as  was  shown  in  works  [1,  3,  4], 
lower  ionosphere  consists  of  three  layers,  C,  D,  and  E.  Layer  C  plays  role  of 
absorbing  pad,  weakening  reflectance  of  layer  D  or  E.  It  was  discovered  as  result 
of  coordination  of  data  on  field  and  medium  by  method  of  mixed  initial  data.  This 
was  done  in  autumn  of  1959  and  reported  at  Colloquium  of  VTs,  Academy  of  Sciences 

of  USSR  on  17  March  i960.  I  take  this  opportunity  to  express  gratefulness  to 
candidate  of  physical  and  mathematical  sciences,  S.  P.  Lomnev,  who  composed 
complicated  program  of  calculations  on  BESM-1  and  carried  them  out  in  short  period. 

7.  Concerning  method  of  normal  waves,  we  note  that  its  application  in  this 
book,  when  wave  numbers  of  normal  waves  are  complex  and  spectra  are  mixed,  would 
have  been  impossible  without  spectral  theory  of  singular  nonself-conjugate  operators 
so  strongly  advanced  during  the  last  few  years  by  Soviet  mathematicians  [16,  J7>  18], 
Now  by  method  of  normal  waves  presented  in  Chapter  II  can  be  solved  problem  of 
waveguide  propagation  of  sound  waves  in  ocean,  infrasonic  waves  in  atmosphere, 
seismic  waves  in  earth,  etc.  Until  now  these  problems,  like  that  considered  in  this 


book,  were  solved  by  quite  complicated  artificial  method  (Watson  method  [7]>  method 
of  mirror  images),  in  which  normal  waves  appeared  without  relationship  to  waveguide 
mechanism  of  wave  propagation.  Moreover,  in  certain  cases  errors  appeared.  For 
instance,  in  work  [19]  were  introduced  waves  THj  and 
and  k.  Besides,  in  the  same  place  and  then  and  in  works  [9,  10,  11]  were  united 
TH0  and  TH^  waves.  In  case  of  application  of  method  of  normal  waves  such  errors 
are  practically  impossible.  We  expended  much  effort  in  order  to  simulate  conditions 
under  which  negative  J  and  k  appear.  For  this  it  turned  out  to  be  necessary  to 
violate  principle  of  radiation  with  respect  to  waves  traveling  on  radius  r. 
Impossibility  of  displacement  of  THq  and  TH^  waves  follows  from  Fig.  10  and  also  is 
due  to  distinction  in  their  forms  in  ideal  case,  when  wallB  of  waveguide  have 
infinite  conductivity. 

Summarizing  what  has  been  said,  one  should  recognize  what  book  holds  interest 
for  specialists  on  wave  propagation  even  at  present.  At  the  same  time  inclusion  in 
it  of  new  ideas  expressed  in  this  preface  would  demand  radical  alternation  of  a 
number  of  its  chapters.  Therefore  it  was  resolved  to  republish  it  by  sterotype, 
using  photomechanical  method. 

25  November  1962 

After  this  preface  was  written,  leading  specialist  on  propagation  of  very  long 
radio  waves,  J.  White,  sent  me  copies  of  his  numerous  works  (see  references  [22-27]), 
illuminating  contemporary  state  of  art  in  the  United  States.  From  these  works  we. 
find  that  starting  from  i960  he  rejected  flat  isotropic  model  of  mrdium  used  in 
works  [20,  21]  and  began  to  consider  sphericity  and  magnetic  field  of  Earth 
approximately  as  this  is  done  in  our  book.  However,  he  applied  Watson's  classical 
method  of  calculation  of  field  [7],  which  is  not  connected  with  contemporary  spectral 
theory  of  linear  operators  [16,  18],  and  results  obtained  by  this  method  pertaining 
to  uniform  routes  (day),  without  explaining  daily  variation  of  very  long  wave  [VLW] 
(CjQB)  field,  in  no  way  correspond  to  results  obtained  in  this  book  or  to  method  of 
calculation  of  fields  (method  of  normal  waves)  used.  Limitedness  of  results  obtained 
in  works  [22-27]  also  is  connected  with  fact  that  their  author  did  not  attempt  to 
coordinate  all  known  data  on  field  and  medium  and  disregarded  factual  information  on 
field  of  -waves,  which  prevented  him  from  constructing  general  picture  of  VLW 
propagation  under  terrestrial  conditions.  All  these  circumstances  once  again  . 
confirmed  my  confidence  in  expediency  of  republication  of  the  book. 


TE^  with  negative  indices  J 
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Author  of  preface  trusts  that  methods  offered  for  coordination  of  data,  method 
of  normal  waves,  and  method  of  modulated  normal  waves,  which  were  used  by  him  for 
formulation  of  theory  VLW  under  terrestrial  conditions,  will  find  application  in 
other  stratified  and  approximately  stratified  media.  He  also  takes  opportunity  to 
express  gratitude  to  N.  A.  Yablochkin,  who  developed  a  number  of  numerical  methods 
for  calculation  of  parameters  of  normal  waves  and  directed  calculation  group  which 
obtained  numerical  results  given  below. 

5  October  19 63 
P.  Ye.  Krasnushkii 
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INTRODUCTION 


CONTEMPORARY  STATE  OF  THEORY  OF  VERY  LONG  WAVE 
VLW  PROPAGATION 

Development  of  theory  [VLW]  (CflB)  propagation  may  be  divided  into  two  periods. 
Basis  of  first  period,  which  began  from  first  stages  of  introduction  of  radio  into 
practice  of  long-distance  communication,  is  diffraction  concept,  advanced  in  1904 
by  Rayleigh  [1 3]. 

In  development  of  theory  of  propagation  of  long  and  very  long  waves  on  the 
basis  of  diffraction  model,  in  which  existence  of  ionosphere  was  disregarded  and 
earth  was  considered  as  uniform  conducting  sphere,  took  part  Poincare'  [14], 

MacDonald  [15],  Nicholson  [16],  Sommerfeld  [17]*  Love  [18],  Marx  [19],  Rybchinskiy 
[20],  and  others. 

Concerning  these  investigations,  in  work  [21],  representing  combined  result  of 
investigations  of  the  first  period,  it  is  stated  that,  "many  years  of  work  of  the 
most  reknown  mathematicians  was  insufficient  to  advance  this  initially  single 
appearing  problem  beyond  the  state  about  which  Nicholson  said  that  of  all 
mathematical  problems  it  alone  has  caused  such  great  divergence  of  opinions."  Not 
until  1918  did  Watson  [22]  succeed,  on  the  basis  of  works  [15]  and  [16],  in  making 
final  conclusion  that  by  the  phenomenon  of  diffraction  of  long  waves  around  clobe 
alone  it  is  impossible  to  explain  fact  of  long-range  propagation  of  these  waves. 

Second  period,  founded  on  waveguide  concept,  began  in  1919  with  work  of  Watson 
[23],  In  this  work  as  model  of  medium  is  considered  uniform  isotropic  semiconducting 
globe,  consisting  in  concentric  uniform  semiconducting  isotropic  ionosphere.  Radio 
waves  are  excited  by  vertical  Hertz  doublet  and  propagate  in  spherical  layer  of 
atmosphere,  located  between  earth  and  ionosphere.  Just  as  microwaves  propagate  in 
contemporary  waveguides. 


-11- 


In  development  of  theory  of  propagation  of  long  and  very  long  waves  on  the 
basis  of  spherical  waveguide  model  with  isotropic  conducting  walls  participated, 
besides  Watson  [23],  Riedback  [24]  and  Preminer  [5],  who  considered  sphericity  of 
earth,  and  also  Kenrick  [25]»  Weirich  [26],  Brekhovskikh  and  Ryazin  [27],  and 
Rudden  [28]  who  considered  it  possible  to  simplify  Watson's  model  by  substituting 
radii  of  earth  and  ionosphere  and  that  conductivity  of  earth  is  infinite. 

Stumbling  block  of  all  above-indicated  works  of  second  period  was  transcendental 
equation  in  cylindrical  functions  for  determination  of  wave  numbers  of  elementary 
waves,  called  normal  waves  by  us,  of  which  is  composed  long-range  field.  During  its 
solution  were  allowed  different  approximations  of  cylindrical  functions,  which  led 
co  considerable  errors.  Another  problem  in  shown  works  was  absence  of  reliable 
data  about  properties  of  lower  layers  of  ionosphere,  which  have  to  be  used  in 
Watson  model  or  its  plane  analog  so  that  theory  yields  simple  answer. 

Certain  authors  managed  to  compensate  defects  of  idealization  and  error  in 
calculations  through  proper  selection  of  parameters  of  model  of  ionosphere.  As  a 
result  of  such  "preparation"  was  created  impression  that  their  theories  were  adequate 
to  reality. 

Example  are  works  founded  on  spherical  model  of  medium,  in  which  is  obtained 
good  agreement  with  Osteen-Kogan  empirical  formula.  Here  the  cited  authors  were 
forced  to  ascribe  lower  boundary  of  ionosphere  of  fantastically  high  conductivity 

-  a  =  1.44.10-"  (COSH  -  1,3.10*  (CCSt), 

-•  -  4.MH*  (COSpj  =  3,6.10*  ( COS» ), 

which  is  at  least  two  orders  of  magnitude  greater  than  real  equivalent  conductivity 
of  very  lowest  layer  of  ionosphere,  responsible  for  long-range  VLW  field  in  daytime. 
In  works  using  plane  waveguide  model  it  was  possible  also  to  obtain  coincidence 
with  separate  experimental  measurements  of  amplitude  of  field  by  means  of  proper 
"preparation"  of  parameters  of  ionosphere. 

In  recent  years  new  works  have  appeared  in  which  is  used  plane  waveguide  model 
of  medium,  taking  into  account  magnetic  field  of  earth  within  the  bounds  quasi¬ 
longitudinal  theory.  Here  was  obtained  coincidence  with  experimental  curve  of  decay 
of  long-range  field  in  spite  of  errors  in  cal  ulation  and  inadmissibility  of  quasi- 
longitudinal  approximation.  In  these  works,  as  in  those  mentioned  above,  two  or 
three  unknown  and  arbitrary  parameters  of  ionosphere  allowed  authors  to  explain,  the 
same  isolated  experimental  fact  through  different  combination  of  numerical  values 
of  these  parameters  and  different  types  of  models.  Such  state  of  affairs  did  not 
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favor  development  of  theory  of  VLW,  and  quite  understandably  there  still  exists 
divergence  of  opinions  on  principal  characteristics  of  waveguide  mechanism  of 
propagation  of  VLW.  For  instance,  there  remain  controversial  the  roles  of  sphericity 
of  earth  and  ionosphere,  of  magnetic  anisotropy  of  ionosphere,  of  erosion  of  lower 
boundary  of  ionosphere,  of  ground  conductivity,  and  also  relief  of  earth's  surface 
and  ionosnhere.  At  the  same  time  the  wealth  of  experimental  material -accumulated 
in  o v$r 50  years  of  practice  of  long-distance  communication  on  VLW  has  turned  out 
to  be  outside  the  domain  of  applicability  of  theory  and  presents  a  quite  disjointed 
picture.  Existing  theories  have  not  been  able  to  explain  even  following  well-known 
and  clear  facts: 

a)  diurnal  variations  in  long-range  VLW  field,  which  are  striking  for  their 
regularity  at  every  fixed  point  of  reception  and  for  their  changeability  as  compared 
to  diurnal  variations  at  different  points  of  reception.  Thus,  for  instance,  there 
remains  the  puzzle  of  4-5-fold  increase  of  signal  strength  at  night  on  Nauen-Tokyo 
line  [1]  and  Hawalan  islands-Moscow  line,  a  two-fold  weakening  of  signal  at  night  on 
Hawaian  islands-Tokyo  line  [1],  a  sharp  15-20-fold  weakening  of  night  signal  on 
Rugby-Moscow  line  (Chapter  V,  Section  J),  etc; 

b)  seasonal  changes  in  strength  of  reception  of  long-range  VLW  stations; 

c)  complicated  dependence  of  field  strength  on  distance  at  night  and  on  routes 
of  varying  illuminance; 

d)  amplification  of  signals  of  long-range  VLW  stations  during  flashes  of 
radiation  of  sun,  accompanied  by  loss  of  shortwave  signals. 
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CHAPTER  I 


FORMULATION  OF  PROBLEM  AND  SELECTION  OF  MODEL  OF  MEDIUM 
§  1.  Mixed  Problem  of  Theory  of  VLW  Propagation 

In  spite  of  impossibility  of  using  existing  waveguide  theories  of  propagation 
of  long  waves  for  purposes  set  in  this  work,  we  subjected  them  tc  detailed  analysis, 
since  it  was  riot  clear  from  above-cited  works  how  and  in  what  measure  offered  models 
of  media  differ  from  reality. 

So  that  analysis  of  models  is  effective  it  is  necessary  to  be  free  from  two 
basic  deficiencies  of  preceding  theories:  1)  arbitrariness  of  selection  of  data  on 
ionosphere  and  2)  inaccuracy  of  algorithm  of  calculation  of  electromagnetic  field 
according  to  data  on  medium. 

The  first  of  these  deficiencies  was  removed  by  newt  formulation  of  problem. 
Usually  solved  is  primal  problem  of  theory  of  waves,  where  electromagnetic  field  is 
determined  in  accordance  with  assigned  properties  of  medium.  In  our  case  properties 
of  medium  are  given  incompletely;  therefore  it  is  necessary  to  solve  so-called  mixed 
problem,  for  which  with  known  and  reliable,  but  not  complete,  data  on  medium  and 
also  from  known  additional  data  on  [VLW]  (CAB)  field  are  determined  unknown  data  on 
medium  and  remaining  VIM  field. 

Since  we  are  interested  ii.  data  on  long-range  VLW  field,  then  as  initial  data 
of  mixed  problem  it  is  expedient  to  take  incomplete  data  about  medium  and  data  about 
short-range  VLW  field.  Result  of  theory  in  this  case  will  be  data  on  long-range 
field  of  VLW  and  unknown  data  on  lower  layers  of  ionosphere. 

Solution  of  mixed  problem  is  carried  out  by  following  method.  All  reliable 
initial  data  on  properties  of  ionosphere,  atmosphere,  and  earth  are  included  in 
model  of  medium,  that  is,  determine  structure  of  function  of  dielectric  constant, 
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which  depends  on  spatial  coordinates,  for  instance,  spherical  r,  6,  and  <p  and  also 
on  several  parameters  a,  £,  7 . a,  v,  z . 

I  PI  I — 1 1«  (M»f :  0»T ....  X,y,z . )  I  !•  ( i .  1 ) 

Permeability  is  considered  equal  to  1  in  CG3.  Fact  of  incompleteness  of  data 
on  medium  is  expressed  in  fact  that  for  part  of  these  parameters  x,  y,  z  ...  their 
numerical  values  are  unknown. 

Let  us  assume  that  also  at  a  number  of  points  of  near  zone 

crA  •>*  ip/)  y«u  3 n 

there  is  known  from  experiment  some  component  of  VLW  field 


_  ,,  ,  -« 


(1.2) 


excited  by  harmonic  currents 

,  7-7|r,l,f),  (1.5) 

4: 

localized  in  near  zone.  Then  by  assigning  I  according  to  (1.5)  and  function  ||f| 
according  to  (1.1),  we  solve  primal  problem  of  theory  of  propagation  of  electro¬ 
magnetic  waves  and  find  component  of  field  E,  as  function  of  T  and  He’ll.  Let  us 
record  this  solution  with  help  of  operator  F 

?M»?) ||7(r',l',f;  —  *,?.**.)  III.  (1.4) 


Since  current  T(r',  9',  cp ' )  in  formula  (1.4)  coincides  with  current  of  formula  (1„3), 
realized  in  conditions  of  experiment,  in  which  were  measured  | E”  j  |  and  (j  =  1,  2, 

3,  .  N),  then  near  field  Eg,  recorded  in  (1.2),  should  satisfy  equation  (1.4). 

Therefore  by  putting  (1.5)  in  (1.4)  we  obtain  2N  equations: 


I  1  ~  ®od FrjAj,y\7 [r'ft'.i')  ;  ||7  (r\r,f;  «,M....xjt,f...)||],  V-5) 

Vjm «rg frtJkpVt p(r'fry);||T(f',rys x#,z...)  ||]  (1.6) 

y=  1,2,3 . n 

for  determination  of  numerical  values  of  parameters  x,  y,  z,  ....  If  number  of 
unknown  parameters  x,  /,  z,  ...  is  equal  to  M,  there  should  be  taken  at  least  M/2 
points  (i'j,  9y  q>j),  in  which  amplitude  |Fj|  and  phase  of  near  field  are  known. 

Solving  system  of  equations  (1.5)  and  (1.6),  we  find  numerical  values  of  x,  y, 
z,  ...  and  putting  them  in  (1.4),  we  determine  from  it  long-range  field  ^ . 

Solution  of  mixed  problem  offered  above  is  depicted  schematically  in  Fig.  1. 
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Solution  of  mixed  problem  ' 

initial  data  of  Calculations  Corollaries 

theory  of  theory 


Fig.  1. 


j 


From  diagram  of  Fig.  1  one  may  see  that  solution  of  mixed  problem  is  replaced 
by  solutions  of  primal  and  inverse  problem,  which  are  based  on  the  same  model  of 
medium,  since  algorithm  of  inverse  problem  is  obtained  from  solution  of  primal 
problem  (this  fact  is  shown  by  dotted  pointer). 

Role  of  primal  and  inverse  problems  in  diagram  of  Fig.  1  is  determined  by 
ratio  of  volumes  of  Initial  data  on  medium  and  field  Eg.  If  all  data  on  medium 
necessary  for  calculation  of  field  E  are  known,  that  is,  unknown  coefficients  x,  y, 
z,  ...  are  absent,  then  no  additional  data  on  field  Eg  need  be  introduced  and 
problem  is  turned  into  primal  one. 

In  other  limiting  case  data  on  medium  are  completely  absent;  then  number  of 
unknown  parameters  x,  y,  z,  ...  becomes  infinitely  large  and  entire  field,  both  near 
and  distant,  should  be  used  as  initial  data  of  problem,  which  degenerates  to  inverse 
problem  of  theory  of  waves. 

In  intermediate  case  (mixed  problem)  number  of  initial  data  on  near  field,  and 
this  means  actual  zone  of  near  field  also  from  which  we  take  these  data,  grows  with 
increase  of  number  of  unknowns  x,  y,  z,  ...,  characterizing  degree  of  our  ignorance 
of  properties  of  medium. 

Basic  formula  (1.4)  was  obtained  by  method  of  normal  waves  and  represents  sum 
of  normal  waves  propagating  on  coordinates  9  and  cp  (Chapter  II).  For  calculation 
of  wave  numbers  of  normal  waves  it  is  possible  to  apply  direct  algebraic  methods 


i 
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which  permit  easy  calculation  of  wave  numbers  with  an  accuracy  of  4-5  digits,  which 
is  necessary  for  determination  of  properties  of  distant  field,  depending  on  difference 
of  wave  numbers  of  first  types  of  normal  waves.  Until  now  wave  numbers  were 
determined  from  transcendental  equation  with  considerably  less  accuracy.  Thus 

method  of  normal  waves  permits  us  to  remove  second  deficiency  of  preceding  theories 

I 

shown  in  beginning  of  this  paragraph.  During  use  of  calculation  diagram  of  Pig.  1 

and  with  exact  algorithms  of  calculation  of  near  and  distant  fields,  all  divergences 

between  corollaries  of  theory  and  experimental  data  on  distant  field  can  be  caused 
only  by  inadequateness  of  model  of  medium. 

§  2.  Analysis  of  Preceding  Models  and  Selection 
ot  rtew  Model  of  Medium 

Method  of  solution  of  problems  of  theory  of  VLW  propagation  presented  in  §  1 
places  in  hands  of  researchers  effective  implement  for  analysis  of  models  of  media. 
Using  this  method  during  1950-1955,  we  analyzed  in  detail  different  models  of  media 
given  in  Table  1.  In  f'rst  vertical  columns  of  this  table  are  enumerated  properties 
of  media  utilized  in  models.  In  upper  line  are  shown  types  of  models  or  authors 

offering  and  using  them.  To  every  model  is  allotted  one  vertical  column  in  which 

plus  signs  denote  considered  properties  and  blanks  or  minus  signs  denote  neglected 
properties  enumerated  in  first  column. 

During  the  study  of  waveguide  models  enumerated  in  Table  1,  we  established 
following  important  internal  regularities  of  models. 

1.  For  calculation  of  near  field  of  VLW  (to  300  km)  it  is  possible  to  apply 
flat  model  of  medium.  Number  of  normal  waves  playing  essential  role  in  near  field 
is  1.5-2  tens;  therefore  representation  of  field  of  near  zone  in  the  form  of  sum  of 
normal  waves  does  not  give  effective  method  of  calculation  of  field.  More 
convenient  in  near  zone  is  beam  method,  which  is  used  in  works  [31,  32,  33,  34  and 
35]  for  determination  of  properties  of  near  field  are  explained  by  interference  of 
one  direct  beam  and  2-3  beams  reflected  from  ionosphere  and  earth.  Beam  method  can 
be  applied  also  at  large  distances  (to  1000  km),  where  it  is  necessary  to  consider 
approximately  Sphericity  of  medium. 

2.  For  calculation  of  distant  field  of  VLW  (more  than  2000  km)  in  range  10-20 
kilocycles  strict  calculation  of  sphericity  of  earth  and  ionosphere  is  necessary,  so 
that  in  distant  field  there  dominate  normal  waves  of  first  numbers,  field  of  which 
is  weakened  near  earth's  surface  from  effect  of  adhesion  of  waves  to  concave  surface 
of  ionosphere  [36,  37],  Small  quantity  of  normal  waves  determining  distant  field 
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makes  expedient  application  of  method  of  normal  waves  for  calculation  of  distant 

field. 

Application  of  plane  models  gives  not  only  quantitative  but  also  qualitative 
divergence  from  reality.  Thus,  for  instance,  cable  wave,  being  basic  in  plane 
model,  is  practically  absent  in  spherical  for  frequencies  f  >  10  kilocycles. 

3.  Disregard  for  conductivity  of  earth's  surface  is  permissible  only  on  routes 
passing  above  marine  surface.  On  long-range  routes  passing  above  land  disregard  for 
conductivity  of  "soil"  can  give  computed  value  of  field  strength  2-3  times  exceeding 
real  value. 

4.  Both  near  and  distant  VLW  fields  are  essentially  influenced  1/  magnetic 
field  of  earth.  Disregard  for  vertical  component  of  magnetic  field  of  earth  leads 
to  considerable  weakening  of  reflection  of  VLW  in  region  of  quasi-Brewster  angles 
as  compared  to  observed  values.  Watson's  model  leads  to  10-20-fold  weakening  of 
distant  field  at  night  as  compared  to  day,  whereas  in  reality  there  is  observed 
2-5-fold  amplification  of  night  field  of  VLW.  In  daytime  allowance  for  influence 

o 

of  magnetic  field  of  earth  is  not  so  essential. 

Quasi -longitudinal  theory  of  Booker  [30],  used  in  work  of  Budden  [29],  is 
inapplicable  for  calculation  of  influence  of  magnetic  field  of  earth  on  distant  VLW 
field,  since  with  gland,  angles  of  incidence  of  waves  hypothesis  of  quasi¬ 
longitudinal  character  of  magnetic  field  is  not  satisfied. 

5.  VLW's  striking  ionosphere  at  large  (distant  zone)  and  small  (near  zone) 
angles  of  incidence  are  reflected  from  ionosphere  at  different  heights.  Therefore 
in  model  of  medium  there  should  be  considered  fact  of  finiteness  of  transition 
region  from  atmosphere  at  different  heights.  Therefore  in  model  of  medium  fact  of 
finiteness  of  transition  region  from  atmosphere  to  ionosphere  should  be  considered. 

6.  For  calculation  of  distant  field  on  routes  of  varying  illuminance  (night  — 
day),  with  varying  character  of  earth's  surface  and  also  under  conditions  of 
changeability  of  ionosphere,  it  is  necessary  to  allow  for  changes  of  height  of 
ionosphere,  its  electrical  properties,  curvature  of  earth,  and  its  conductivity  and 
relief  along  route  of  radio  waves  connecting  corresponding  points. 

Our  problem  consists  in  formulation  of  theory  of  VLW  propagation  which  ade¬ 
quately  explains  on  the  basis  of  single  model  of  medium  regular  phenomena  of  near 
and  distant  fields  under  conditions  of  propagation  on  routes  of  mixed  illuminance 
(night  -  day),  varying  properties  of  earth's  surface  (land  —  sea),  and  also  random 
phenomena  caused  by  flashes  on  sun,  by  magnetic  storms,  and  stationary  turbulence  of 
ionosphere  at  night. 
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Considering  in  light  of  these  requirements  imposed  on  theory  of  VLW  the  six 
above-enumerated  points  characterizing  VLW  propagation  in  different  models.  Table  1, 
we  are  convinced  that  in  view  of  nonfulfillment  of  points  4,  5  and  6  models  of  media 
applied  earlier,  in  particular  the  Watson  and  Budden  models,  cannot  be  used  for 
formulation  of  adequate  theory  of  VLW.  In  connection  with  this  we  considered  new 
models  of  media,  designated  in  Table  1  by  letters  A,  B,  and  C. 

Model  A  is  generalization  of  Watson's  waveguide  model  for  case  of  calculation 
of  vertical  component  of  magnetic  field  of  earth. 

Model  B  differs  from  model  A  in  that  in  it  there  is  considered  erosion  of  layer 
of  ionosphere,  which  is  represented  in  the  form  of  series  of  concentric  layers  with 
radii  rk,  k  =  1,  2,  ...  N,  each  of  which  is  characterized  by  tensor  of  dielectric 
constant  ||e.  ||,  remaining  constant  within  limits  of  layer.  Model  A  is  particular 
case  of  model  B  for  N  -  1. 

Model  C  represents  further  generalization  in  this  series  of  models.  In  it  is 
considered  fact  of  slow  changes  of  radii  of  layers  r^  and  their  electrical  properties 
117^11,  depending  upon  geographic  coordinates  B  and  <p. 

From  6  above  mentioned  points  it  follows  that  only  model  C  satisfies  all 
requirements  imposed  on  VLW  theory.  Model  B  has  field  of  application  (uniform 
route  under  night  or  day  conditions) .  Model  A  can  be  applied  only  for  approximate 
calculations,  since  it  does  not  permit  us  to  use  method  of  introduction  of  initial 
data  of  near  field  as  per  diagram  depicted  in  Fig.  1. 

It  is  necessary  to  note  that  in  region  of  near  field,  which  we  shall  use  in 
Chapter  III  during  introduction  of  initial  data,  model  C  is  replaced  by  model  B. 

Thus  as  initial  data  on  Eg  are  introduced  averaged  quantities;  therefore  unknown 
parameters  of  medium  and  He^ll  obtained  from  theory  also  are  averaged  over 
distances  of  300-500  km. 

§  3.  Plan  of  Work 

As  basis  of  this  work  are  models  of  media  B  and  C  described  in  preceding 
paragraph.  Solution  of  mixed  problem  of  theory  of  propagation  of  waves  is  carried 
out  by  method  described  in  §  1.  Therefore  diagram  depicted  in  Fig.  1  is  essentially 
plan  of  present  work. 

In  Chapter  II  per  diagram  of  Fig.  1  are  constructed  models  of  media,  anu  by 
method  of  normal  [36,  37)  waves  is  solved  primal  problem  of  theory' of  propagation 
of  radio  waves  for  models  B  and  C. 
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In  Chapter  III  are  given  results  of  solution  of  inverse  problem  of  theory  of 
waves,  determining  unknown  parameters  of  lower  layers  of  ionosphere  x,  y,  z,  .... 

In  Chapters  IV-VI  we  return  to  primal  problem,  introducing  in  it  known  data  on 
medium  (a,  0,  y,  . ..)  and  new  data  on  ionosphere  (x,  y,  z,  ...)  obtained  in  Chapter 
III. 

Chapter  IV  contains  basic  characteristics  of  normal  waves  and  their  dependence 
on  parameters  of  medium  and  radiator. 

Chapter  V  is  devoted  to  calculation  of  regular  processes  of  VLW  propagation 


over  large  distances 


CHAPTER  II 


SOLUTION  OF  BOUNDARY  PROBLEM  OF  VLW  PROPAGATION 
ON  BASIS  OF  MODELS  B  AND  C 

§  1.  Introduction  of  Data  on  Properties  of  Medlm 
Represented  by  Models  B  and  C 

Models  B  and  C,  briefly  described  in  §  2,  Chapter  I,  have  so  general  form  that 
they  can  be  applied  to  any  radio  waves  propagating  around  the  globe  by  means  of 
ionosphere  when  it  is  possible  to  disregard  small-scale  turbulences  and  horizontal 
component  of  magnetic  field  of  earth  (high  and  middle  latitudes).  If  we  apply  models 
B  and  C  in  so  general  form  for  solution  of  mixed  problem  per  diagram  of  Fig.  1,  large 
quantity  of  experimental  data  on  field  E  (near  and  distant)  will  be  demanded,  which 
have  to  be  introduced  as  initial  data  of  theory,  so  that  theory  yields  simple  answers. 
Here  the  main  part  of  corollaries  of  theory  will  consist  of  information  represented 
by  models  B  and  C,  and  only  small  part  will  represent  distant  field  E.  Considering 
that  our  main  aim  is  to  obtain  new  information  about  distant  field  E  and  not  study 
of  properties  of  lower  layers  of  ionosphere,  we  should  make  models  B  and  C  as 
specific  as  possible  by  using  sources  of  information  unbound  to  experimental  data 
on  [VLW]  (c,QB)  field.  Here  problem  will  approach  the  primal,  and  volume  of  additional 
experimental  information  about  field  Eg,  utilized  as  initial  data  of  theory,  will 
decrease.  We  have  all  data  necessarv  to  represent  two  first  layers  (0,  rQ)  and 
(*0*  *!.)• 

Earth  layer  occupies  interval  (0,  Vq) ,  where  rQ  =  a  is  radius  of  earth's  surface. 
Jince  VLW's  do  not  penetrate  earth  or  water  by  more  than  1-2  km,  we  shall  consider 
that  whole  globe  possesses  uniform  conductivity  a  =  Oq  and  is  characterized  for 
frequency  u>  by  complex  dielectric  constant. 

<>  •  »  (2*1) 
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corresponding  to  effective  value  e0  in  upper  layers  of  earth  for  frequency  f  = 

=  ^  =  10-20  kilocycles.  Numerical  values  of  rQ  and  Eq  are  well  known  and  will  be 
introduced  later,  in  Chapters  IV-V. 

Following  layer  of  models  B  and  C  is  atmosphere  (Tq,  r^) .  From  data  on 
measurement  of  dielectric  constant  in  lower  layers  of  atmosphere  it  follows  that  for 
VLW  dielectric  constant  of  layer  of  atmosphere  may  be  considered  equal  to  i. 

Data  on  layers  of  ionosphere  (i^,  r?) ;  (r2,  r-j) ;  (r-j,  r^) . (’qj-l*  ®)  are 

scanty  and  little  reliable.  From  literature  known  to  us  about  lower  layers  of 
ionosphere  we  are  unable  to  obtain  reliable  quantitative  data  about  radii  and 
tensors  ||ejJI  (k  =  2,  3,  ...  N) .  Together  with  this,  from  these  sources  it  is  possible 
to  represent  structure  of  ionospheric  layers,  which  considerably  decreases  number  of 
necessary  initial  data  of  theory  of  field  . 

The  most  reliable  pieces  of  information  about  lower  layers  of  ionosphere  are 
the  following. 

1.  Experimental  data  on  absorption  and  reflection  of  long,  medium,  and  short 
waves,  indicating  existence  in  daytime  of  ionospheric  layers  below  E  layer  and 
absence  of  layers  of  ionization  below  E  layer  at  night  [1-5]. 

2.  Theoretical  research  [38,  39]  on  microprocesses  in  lower  layers  of  ionosphere, 
from  which  it  follows  that  at  heights  of  70-90  km  basic  role  on  frequencies  of  10-20 
kilocycles  is  played  by  electrons. 

3.  Theoretical  treatment  of  different  experiments  [40],  allowing  us  to  estimate 
dependence  of  effective  number  of  collisions  veff  electrons  with  other  particles  on 
height  of  layers  h  in  interval  of  60-100  km  (Fig.  2) . 

4.  Experimental  data  on  magnitude  of  magnetic  field  of  earth  HQ. 

From  first  point  above  the  conclusion  can  be  made  that  very  long  waves  are 
reflected  from  ionospheric  layers  located  at  heights  of  70-90  kir.  which  it  is  possible 
to  approximate  by  concentric  (model  B)  or  approximately  concentric  (model  C)  layers 
(r1,  r2),  (r2,  r?)  ...  (rN_1,  co),  where  (rx  -  rQ)  =  70  km  and  (rN_1  -  rQ)  =  100  km. 

Thus  the  last  layer  of  model  of  ionosphere  (r^^,  °°)  io  propagated  infinitum. 
Distinction  between  model  of  medium  from  real  ionosphere  in  region  r  >  ^  is 

immaterial,  since  VLW's  do  not  penetrate  of  E  layer. 

Another  important  point  for  accurate  representation  of  model  of  medium  is 
electron  character  of  conductivity  of  luwer  lasers  of  ionosphere.  Considering  that 
motion  of  electrons  under  action  of  harmonic  field  Ee"*^  in  lower  layer  of  ionosphere 
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is  described  by  Lorentz  equation 

—  immV  + 

C 


where  to  —  angular  frequency  of  electromagnetic  field, 
m  —  mass  and  e  -  charge  of  electron, 
v  —  velocity  vector  of  electron, 

we  obtain,  disregarding  horizontal  component  Hq,  following  tensor  of  dielectric 
constant  for  k  layers 


for 

where 


I*  °*  °  | 

°,  % 

fi-i  (A  ■  2,  3,  4  N), 


•ff- 


1  + 


_ H-<y _ 

•Mv-HF+«W 


(2.2) 

(2.5) 

(2.4) 

(2.!?) 


In  formulas  (2.3) -(2.5)  are  introduced  following  conventional  designations: 


(2.6) 

M 

mH - - gyromagnetic  frequency. 

(2.7) 

Index  of  k  layer  in  formulas  (2. 3) -(2. 7)  is  omitted.  Thus  every  "k"  layer  of 

? 

ionosphere  is  characterized  by  four  quantities  o)Q,k,  u^,k,  and  r^. 

Quantity  vg^,k  is  connected  with  height  of  ionosphere  h  =  -  rQ,  according 

to  work  Nlcolet  [40],  by  formulas 


which  are  depicted  in  Fig. 


v"m*ipFF  for  A  <90  km, 

Vf  =  V4.1(p?p  for  A  >  80  km, 

2. 


-25- 


A*- 


Quantity  is  calculated  from 
experimental  value  of  and  is  approximately 
equal  to  8.10^  for  all  layers  in  interval 
of  70-90  km  at  mid  and  northern  latitudes. 
Thus  unknown  quantities  in  models  B  and  C 
are  distribution  of  electron  concentration, 

p 

determined  by  quantities  u>q  ^  and  r  . 

Further  specification  of  models  B  and  C  on 
the  basis  of  known  data  on  medium  is  not 

p 

possible,  and  parameters  cack  and  r^,  which 
we  designated  in  Chapter  I  through  x,  y, 
z,  ...,  are  determined  from  solution  of 
mixed  problem  per  diagram  of  Fig.  1. 


§  2.  Formulation. of  Boundary  Problem  on  Basis  of  Model  B 

In  accordance  with  method  of  solution  of  mixed  problem  presented  in  §  1, 

Chapter  I  we  first  solve  primal  problem  of  theory  of  VLW  propagation  on  the  basis  of 
model  B,  considering  that  all  parameters  of  model  B  are  known.'  Primal  problem  of 
VT.W  propagation  on  the  basis  of  model  B  is  formulated  in  the  following  way. 

In  medium  there  exist  distributed  outside  currents,  harmonic  in  time  and  directed 
radially 

t  (2.8) 


It  is  assumed  that  distribution  function  of  outside  currents  Ir(r,  0 )  depends 
only  on  spherical  coordinates  r  and  9,  that  is,  is  characterized  by  symmetry  of 
rotation  with  respect  to  polar  axis  0  =  0  of  spherical  system  of  coordinates  r,  6,  9. 

To  be  found  is  stationary  electromagnetic  field  created  by  currents  (2.8)  in 
spherically  laminar  medium  described  by  following  tensors  of  dielectric  constant 


IJhH  = 


0 

0 


.0,0 

.•A*!? 


(2.9) 


where 


for  zero  layer  (k 


r^t<r<rt\  *  =  0, 1,2,3... Af, 

0);  {r_1  ---  0,  r0) 


(2.10) 


for  first  layer  (k  =  1);  (rQ  =  a,  r^) 


•rr  — «it  — 0;  (2.11) 

and  for  other  N  - 1  layers  (k  =  2,  3,  ...  N) 

(r*»  r$)i  (rtt  **)•••  [rN-i>  oo)  (2.12) 

tensor  He^ll  is  described  by  formulas  (2.3)-(2.7). 

In  virtue  of  symmetry  of  excitation  currents  (2.8),  electromagnetic  field  should 
be  determined  from  Maxwell  equations,  in  which  all  derivatives  with  respect  to  are 
equal  to  zero.  Thus  for  every  spherical  "k"  layer  components  of  field  E  ,  E0,  Er, 
and  H^,  H0,  Hr  are  determined  by  following  system  of  Maxwell  equations,  recorded  in 
spherical  system  of  coordinates  (CGS): 


(2.13) 

ry  >rW*)-  t!1**** + *•»  *!*• 

(2.14) 

-  “T  +  w4*1 ' 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

in  interval  rk-1  s  r  s  k  =  0,  1,  2,  3,  ...  N. 

On  boundaries  of  layers  have  to  be  met  conditions  of  continuity  of  tangential 
components: 


(2.19) 

(2.20) 

=«•/%. 

(2.21) 

(2.22) 

where  k  =  0,  1,  2,  ...  N  -  1. 

Besides,  there  have  to  be  met  conditions  of  limitedness  of  solutions  at  point 
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/* 


| 
f  ' 


and  condition  at  infinity 


mod  £|/-,  <  At;  mod  HJr_,  <  M 


mod  ZJr -*•.-*■  0;  mod  /7}r-*»  -♦0. 


(2.23) 


(2.24) 


Equations  (2. 13) -(2. 18)  under  conditions  on  boundaries  of  layers  (2,19) -(2.22) 
and  conditions  at  singular  points  (2.23) ,  (2.24)  for  finite  values  vef,f  N  and  Ng  ^ 
in  simple  form  determine  electromagnetic  field  (E,  H)  for  given  T 


crop  * 


§  3.  Introduction  of  Electrical  and  Magnetic 
Potentials  A  and  B  of  Potential 


Vector  Function 


For  solution  of  boundary  problem  formulated  at  the  end  of  preceding  paragraph 
we  introduce  two  scalar  functions  A  and  B,  determined  through  components  and  H^: 


_  1  dA 

r'-f  ^  —  > 

r  dA 

mj  I  dB 

"’"rfi- 


(2.25) 

(2.26) 


Then,  according  to  ( 2. 13) -(2. 18) ,  all  remaining  components  are  determined  through  A 
and  B  in  the  following  way: 


„  l  d  f ,  ,dB\ 

jk,«rrr* slnO  dH  \  "  d8/’ 


iPB 


*0? 


dA 


Ai»66r  drd'i  r<Vl  d O’ 


i 


A,  r*  sin  8  <W 


l-S) 


«  1  *A 
H%~ktrdrd*' 


(2.27) 

(2.28) 

(2.29) 

(2.30) 


where  K1  =  —  is  wave  number  of  free  space. 

After  introduction  of  functions  A  and  E,  system  of  Maxwell  equations  (2. 13) -(2. 18) 


4 


takes  f“rm  of  series  of  N  +  1  pairs  of  connected  scalar  wave  equations  of  second 
order: 


r +  oti  n  (i,n  0  [' 1# +  A' 


Iff 


(2.31) 

(2.32) 
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v  A 


in  interval  .  Sr"  r^  k  -  0,  .,  2,  ■ ,  ...  N. 

System  ol'  equations  (2.31)  and  (2.  ),  according  to  (2. 19) -( 2. 22) ,  should 

satisfy  conditions  of  jumps  for  A  and  1*  on  internal  boundaries  of  media 


(2.5?) 


(2.34) 


(2.5?) 


(2.?6) 


where  k  =  0,  1,  2,  3,  .  .  .  N  -  1 , 

and  also,  according  to  (2.23)  and  (2.24),  boundary  conditions! 

modA9<M'\  mod  B%<N'  for  r-*0, 


(2.37) 


where  M'  and  N'  are  finite  real  numbers 

otodi4-»0;  mod B -*■ 0  f n>  r -*•  oo. 

Now  we  introduce  vector  function  of  potentials  with  components  A  and  B 


vector  function  of  currents 


(2.38) 

and 


.  Then  we  record  system  of  equations  (2. 31) -(?. 52) 


in  operator  differential-matrix  form 

Bu 

A* 


♦mi-Ma. 


(2.39) 


where  is  matrix  operator 


/<*>  = 
t 


'*•%  i 


dr  * 


*M 


d 1 


it  *■ 


[<*¥} 


(2.40) 


in  interval  (r^.^.  r^)  k  =  0,  1,  2,  ...N, 
and  l^k)  i-  matrix  operator 


-2d- 


IKPIIIKJWI*!  Wt*  - 


<!-/, 


1  *1 

(sin  4  5T 


Corresponding  conditions  of  jump  take  form: 

1  d 


IT 


± 

dr 


*>) 

0 

1 

tin 

take 

form: 

1  d 

* 

ftS 

>'* 

0 

Ifl.l 

-M 

lv„ 

s*s(,l"'s)- 


(2.41) 


**,•£►*» 


‘IT 

L- 

dr 


B 


*+i 

U 


(2.42) 


(2.43) 


*-=0,1,2  ....  N-\, 

and  boundary  conditions  (2.37)  and  (2.38): 


mod 


CM? 


wh»n  r  -*•  0, 


U 

mod  “  >0 


-htrr 


or.. 


(2.44) 

(2.4*,) 


Since  differential-matrix  operators  determined  by  (2.40)  differ  from  each 

/  lr  \  lr  Ir 

'then  only  in  numerical  values  of  parameters  e^,r  »  eqq>  and  eG cp’  wa  consider 

system  of  equations  (2.39)  k  =  0,  1,  2,  ...  N  for  intervals  (r^_j_»  r^),  k  =  0,  1, 

2,  ...  N  as  one  equation  in  interval  (0,  ai) 


A  *  A  c  0 


(2.46) 


in  which  operator  l  consecutively  takes  forms  .  ...  in 

r  r  r  i*  r 

accordance  with  passage  through  varying  r  intervals  (0,  Tq) ,  (r0,  r^),  (r^,  r2) 


B* 


,  k  =  0, 


(ri,  r2)  ...  (rN_1#  co ).  Vector  function  |ft|,  taking  values 

1,  2,  ...  N,  in  corresponding  intervals  (0,  r0),  (rQ,  r^)  ...  (rN_1,  oo)  undergoes 
breaks  at  points  rQ,  r1,  r2,  ...  rN_1,  where  it  satisfies  conditions  of  break 
(2.42)-(2.43). 


For  N  +  i  layers  (0,  rQ) ,  (rQ,  Tj),  ...  (r^,  ®)  there  are  2N  conditions  of 
break  (2.42)  and  (2.43)  and  two  conditions  (2.44)  and  (2.45)  at  boundaries  of  interval 
(0,  co),  2(N  +  1)  conditions  in  all,  determining  2(11  +  1)  arbitrary  constant  solutions 
of  equation  (2.46). 
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Operator  10  has  two  singular  points  0=0  and  9  =  v  at  which  we  shall  require 
limitedness  of  solution  |?|  of  equation  (2.46). 


&■ 

t 


§  4.  Operator  of  Normal  Waves  and  Its  Spectrum. 

Free  formal  Waves 

For  finding  solution  of  nonuniform  equation  (2.46)  under  conditions  (2.42) -(2.46) 
we  apply  method  of  normal  waves  presented  in  works  [36,  37].  According  to  method  of 
normal  waves,  we  first  find  particular  solutions  of  homogeneous  equation  (2.46)  for 
conditions  (2.42) -(2.45) ,  having  form 

/  |fl|+  <|  flL0. 

'U I  iU  (2. '17) 

Let  us  find  particular  solution  of  (2.47)  in  the  form 


nr) 

Z[r ) 


*(•). 


where 


W  _  K,(r) 
Z(r)  Zk[r) 


in  interval 

r*-i  <r<rte 

Putting  (2.48)  in  (2.47)  and  dividing  variables,  we  obtain 

nn  | 


(2.48) 

(2.49) 


l. 


Hr)  |  _ 

Z[r)  | 


(2.90) 


Identity  (2.50)  can  exist  only  under  the  condition  that  both  its  parts  are  equal 
to  the  same  constant,  which  we  designate  by  letter  Then  we  obtain  two  equations: 

™Uzl 


nr) 
Z(r)  r 


Z(r) 

/0r(0)  +  XC(0)~0 


(2.51) 

(2.52) 


First  equation  in  expanded  form  will  be 


Y 

Y 

X 

Z 

-Z* 

Z 

-*trV  +*V»+ 

(2.55) 


| 

! 

1, 


i 

* 

I 


I 

I 

X 


I 

f. 


where  err  =  e*r;  eQ6  =  e£0J  e0q)  =  e*,  in  intervals  (r^,  rR) 

for  k  =  0,  1,  2,  3,  . . .  N. 
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On  boundaries  of  interval: 


V  k 


0,  1,  2,  3,  ...  N  -  1  conditions  of  Jump  are  met: 


1  d 

0 

•w 

d_ 

X 

Y. 

Zu 

dr 

i 

1 

1  d 

| 

- 

0 

•IV11  * 

d 

IF' 

X 

Yu+i 
Zu*  1 

,th 


Y> 

Z> 


Yu  f. 
Zu+i 


(2.54) 

(2.55) 


and  on  ends  of  zero  and  N  intervals  boundary  conditions 


mod 

mod 


HL- 


limit 


(2.56) 

(2.57) 


Differential-matrix  expression  presented  in'  (2.53)  for  boundary  conditions 
(2.54),  (2.55/1  (2.56),  and  (2.57)  is  operator  of  normal  waves,  which  acts  on 
vector  function  [Y/Z].  Designating  it  by  symbol  Lr,  we  record  equation  (2.53)  with 
conditions  (2.54) -(2. 57)  in  the  form 


\z\ 


(2.58) 


Equation  (2.58)  determines  spectrum  of  operator  of  normal  waves.  In  general 
operator  Lr  is  not  self-conjugate  and  singular.  Its  spectrum  consists  of  discrete 
and  continuous  parts.  In  case  of  great  conductivity  of  earth  |f.q|  »  1,  solid  part 
of  spectrum  may  be  disregarded  in  examining  field  in  interval  (rQ,  id)  and  under  the 
condition  that  exciting  currents  are  located  outside  earth.  Therefore  we  shall 
consider  subsequently  only  discrete  part  of  spectrum,  that  is,  eigenvalues  of  operator 
Lt,.  Let  us  renumber  them  in  order  of  growth  of  absolute  values 

*•»  *!»  X.  (2-59) 

Since  operator  L  is  not  self-conjugate,  eigenvalues  of  will  be  complex 

t, -R*(  X,)  +  Um[  Xj  )  y-0,1,2...  (2.59') 


From  what  is  to  follow  ne  will  sec  that  to  every  eigenvalue  of  Xj  there 


corresponds  only  one  eigenvector  function 


ea-1^ 


for  rk-1  s  r  s  rk  k  =  0,  1,  2,  ...  N. 


(2.6o) 


Eigenvector  functions  (2.60)  of  operator  Lr  are  orthogonal  eigenfunctions  of 

m 

(4¥HM). 


conjugate  operator  Lr*,  which  is  determined  from  Lagrange  identity 


where  parentheses  designate  scalar  product  of  two  vector  functions. 

Let  us  designate  eigenvalues  and  vector  functions  of  conjugate  operator 
respectively 

N»  Ihi  l*i  •••  l1^  ••• 


which  satisfy  equation 


tm  •••  n  - 

From  (2.58)  and  (2.65)  under  condition  (2.6l)  follows  relationship 

IS  I  MU 0  for  y*. 


(2.61) 


(2.62) 


(2.65) 


(2.64) 


or  in  scanned  form 


00 

J" 


(r)(/,(r)rfr  +  j*Z/(rj.  V,{r)dr-0;  j  +  p , 


(2.64*) 


expressing  orthogonality  of  eigenfunctions  of  conjugate  operators  Lr  and  Lr*. 

For  operator  I.  ,  generated  by  expression  (2.55)  and  boundary  conditions 
(2.54) -(2.57) ,  it  is  easy  to  show  that 


(2.65) 

(2.66) 


where  vinculum  denotes  conjugate  complex  numbers. 

In  virtue  of  (2.65)  and  (2.66),  condition  of  orthogonality  for  non-self -conjugate 


operator  Lr  takes  form 


BIBFD- 


(2.67) 
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or  in  expanded  form 


-o 


If  j  =  P»  then  (2.64)  is  equal  to  Nj,  called  subsequently  the  normalising  factor 

+  (2.69) 


Now  let  us  return  to  second  equation  (2.52),  which  for  eigenvalue  of  Xt  in 


expanded  form  has  the  form 


Introducing  new  designation 


(2.52') 


(2.70) 


we  convert  it  to  Legendre's  equation,  which  has  two  particular  solutions  [37]: 


(2.71) 


(2.72) 


:s  particular  solution  (2.48)  of  equation  corresponding  to  eigenvalue  of  Xj 


ha:-  the  form 


(2.7-i) 


Here  we  introduced  time  factor  e  .  In  region  of  permissibility  of  asymptotic 
^presentations  Lf 1 ^  and  [37]  we  record  (2.73)  in  the  form 


(2.7*0 


where 


(2.75) 


(2.74)  represents  suiu  of  two  waves  traveling  in  positive  and  nera'ive  directions  on 
coordinate  ■1  with  angular  phase  velocity  determined  by  angular  wave  number  v  ..  both 
waves  are  modulated  in  identical  manner  with  respect  to  front  of  wave  (•  -  const 


according  to  eigenfunction  |Yj( r) | Zj(  r)  |  .  Such  waves  were  called  normal,  since  they 
were  obtained  as  a  result  of  solution  of  problem  for  eigenvalues  of  operator  of 
normal  waves,  just  as  normal  oscillations  are  obtained  as  a  result  of  solution  of 
problem  for  eigenvalues  of  operator  of  normal  oscillations. 

Sum  of  normal  waves  (2.75)  represents  any  wave  process  harmonic  in  time  in 
region  where  exciting  currents  3^  are  absent.  Therefore  by  analogy  with  normal 
oscillations,  such  waves  are  called  free  normal  waves  [56,  57]. 

Solution  of  equation  (2.46)  on  right  side  yields  forced  normal  waves,  which  v;i 
be  considered  in  the  following  paragraph. 

§  5.  Formulation  of  Solution  of  Nonuniform  Boundary 
Problem.  Forced  Normal  Waves 


Now,  following  method  of  normal  waves,  let  us  look  for  solution  of  nonuniform 
equation  (2.46) 


L, 


B 


(2.76) 


in  the  form  of  spectrum  of  normal  waves 


B 

A 


2 


YA') 

ZAr) 


(2.77) 


We  recall  that  Lr  is  operator  with  respect  to  coordinate  r,  generated  by 


differential  matrix  expression  (2.55)  and  boundary  conditions  (2.54)-(2.57),  while 
L q  is  operator  with  respect  coordinate  0,  generated  by  differential  expression  (2.41) 
and  conditions  of  limitedness  of  solution  at  points  0  =  0  and  0  =  tt. 

For  determination  of  unknown  phase  factors  we  put  (2.77)  in  (2.76) 


(2.78) 

\Z,  I  "  IZy  I  C  10  1 

M  /-« 

Sine  series  (2.77)  converge  evenly,  we  change  order  of  summation  and  operations  Lr 
and  Lr 


"0 


y‘  +  y  r‘  *,  */--'* 

TZ  zj  ri  Zj  *  c 


(2.79) 


Multiplying  both  parts  of  (2.79)  scalarly  on  some  eigenvector  of  conjugate  operator 


Lr*  and  taking  into  account  (2.58),  we  obtain 


(2.80) 
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In  virtue  of  ronditi  n  of  orthogonality  of  (  '.<4),  expression  (2.80)  is  .  imp !  i  •  ltd 


(  .  t"V 


where  N  is  normalizing  factor  (2.69). 

Expanding  operation  L(  and  using  (2.66),  we  obtain  equation 


1  L( 

•Int  Jl\ 


tint 


*•< 

~dV 


9 


+*j*n 


(>v  ) 


which,  under  condition  of  limitedness  of  solution  at  points  0  =  0  and  0  n , 
determines  phase  factors  $j(0)  °t'  normal  waves.  These  factors  differ  from  phase 
factors  ¥j(0)  of  free  normal  waves,  considered  in  §  4,  Chapter  II,  since  they  an 

forced  solutions  of  Legendre  equation  on  right  side  (2.81). 

F’orced  solution  of  equation  (2.81)  can  be  found  by  different  methods.  In  work 
1^7]  it  was  obtained  with  help  of  Green's  function.  Here  we  shall  solve  equation 
(2.81)  by  method  of  variation  of  arbitrary  constants  of  general  solution  of  homo¬ 
geneous  Legendre  equation,  which  it  is  possible  to  represent  in  the  form 


A,PV  f  coi(«  -  #)j  +  BjPv(  cos  6). 


Considering  and  B,  functions  of  0,  we  find  them  from  two  following  equal  i 


*/<wy  -« +«/«)*,(*> =0, 

vwy-l)  +  */(l|PvI«l  -  ~rl,- 


where  primes  denote  derivatives  with  respect  to  0,  and 

I— cost, 

W* 


_  f  AW 


From  these  equations  it  follows  that: 


4,  1 

*•»-- tstJ— i 5— 


2  tin 
«  ilnl 


(".>  '8 
C  •  '-'1 1 

( -V  •  1 
(■ .  * ) 


(  . 


where 


Thus  particular  solution  of  nonuniform  Legendre  equation  (2.81)  has  the  form 

..  4.7  t'Af-W 


and  general  solution,  after  substitution  of  (2.88), 


•^«)-C/>,|co.(.  -  4)1  +  D,  P„|co.4|  - 

*,  i 

+  ~nV>v.  M^l 1  j '/(*')/\;|co.(.-4')|.mW . 


f! 

V 


It  should  satisfy  conditions  of  limitedness  of  function  $j(£)  at  points  0=0 
and  9  =  rr.  Fulfillment  of  this  condition  at  0  =  0  occurs  only  if  Cj  =  0.  Indeed, 
when  0=0 

•/(I)  -  CjP^[  coi «)  +  D/^coi  0). 

Since  function  Pvj[cos(7r  -  0)]  with  0  -►  0,  has  logarithmic  peculiarity,  it  is 
necessary  that  Cj  =  0. 

At  point  0  =  tt  general  solution  will  have  form 


« 

— A/co,0)J/;(4-)A;|«.>-)«lnlW. 


•/(*)— A(C0««) 


It  will  be  limited  at  this  point  if  coefficient  for  Pv(cos  ir)  is  equal  to  0, 
whence 


ft- 


Thus  solution  of  equation  (2.81)  will  be 


'A  i«r  )•*«' 


(2.89) 


% 

fi 


* 

f 

1 


I 
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r 


(2.09 
cont 'd) 


Complete  solution  of  initial  equation  (2.76),  according  to  (2.77) ,  will  have 


form 


W I _ L 


B  --—Vi 

A  e  ijzy(r)|  ilnyt-A^ 


/yco»(«-i)|x 


X I P,ytco.  I>m  IW  +  /y  cos  IJ  •  J  //yco*  -  r))iln  W 


(2.90) 


In  particular  case,  when  field 
on  polar  axis  0  =  0  at  point  r  =  b. 


B| 


is  excited  by  vertical  Hertz  doublet,  located 


|»| <d|  >/W-.p,[cOT(«-.)|, 

Ml  rt*  “Iz.wl  / 


(2.91) 


where  P  is  electrical  moment  of  dipole. 

Phase  function  0j(0)  In  this  case  has  the  form 

#  = 


►,= -■  Jv', 

•  #*cA(y*ln»y«  7 


(2.89*) 


In  order  to  consider  field  at  large  distances  0  from  Hertz  doublet,  we  use  formulas 
(2.71)  and  (2.72)  and  present  P^,[/-os(fr  -0)]  in  the  form 

*»,!«•(«  -  •)  I  =  ~  { L^lcoK*  -  •)  J  -  l)jj. 

i 

Further  we  expand  ~  in  series  after  formula  of  geometric  progression 

J 


*■*,« 


=  -2 It  1  7  +  ..•* 


jl  +  ...J. 


2irv  .i 


It.  I'  K  *  X  j 

e  J  <  1,  if  we  take  principal  values  of  roots  Vy 
Thus  we  present  (2.91)  in  the  form 


Switching  to  asymptotic  representation  of  functions  j'  and  L^y,  which  is 


Justified  for  sufficiently  large  0,  we  obtain 

-U[B\  2 P 


J£  .  /  n'Vf/A 

^5*  J/  iitl  W|z,W|T87*  X 
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(2.92) 


From  (2.Q2)  one  may  see  that  second  members  in  brackets  represent  normal  waves 
arriving  at  point  of  reception  (r,  0)  from  radiator  0=0.  Where  wave  n  =  0  passes 
over  the  shortest  arc  of  great  circle,  lagging  in  phase  by  Rg(Vj*)  0,  and  waves 
c  n  i  0  pass  n  times  around  globe  reaching  point  (r,  0);  they  are  called  direct 
echoes. 


First  members  in  brackets  are  normal  waves  coming  to  point  of  reception  (r^  0) 
over  larger  arc  cf  great  circle  (n  =  0),  having  first  circled  the  globe  n  times 
(inverse  echoes). 

If  Vj  have  sufficiently  large  imaginary  parts,  direct  and  inverse  echoes,  and 
also  first  reverse  wave,  are  considerably  weaker  than  normal  waves  coming  to  point 
uf  reception  by  the  shortest  means,  and  it  is  possible  to  disregard  them.  Then 
approximate  solution  (2.92)  we  record  in  the  form 


B 

A 


(2.93) 


§  6.  Calculation  Formulas  for  Components  of  Electromagnetic 
Field  in  Layer  of  Atmosphere  and  Under  Surface 
of ' Earth  (Model  b) 

We  shall  use  (2.93)  subsequently  for  calculation  of  electromagnetic  field  in 
layer  of  atmosphere  (rQ,  r^)  and  under  water  or  earth's  surface  (r  <  rQ) .  In  both 
these  layers  media  are  Isotropic;  therefore  components  of  electromagnetic  field  are 
calculated  by  simplified  formulas  (2.29) -(2.30) ,  with  =  -e^g  =  0. 

For  every  normal  wave  of  number  J 


r  .*  “  r  (ft 


r  (ft 


*,r* 


(2.9*) 
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e  l  9B,_  Idyjd+i.u  -  1  d'A/  _  1  dZJ  d*t  (2.94 

•/“  A,r  drdt  ~  kxr  dr  '  dt  ’  •/  kf  drdb  *,r  dr  '  <rt  '  cont'd) 

Here  Is  used  equation  (2.52). 

As  can  be  seen,  components  of  field  are  divided  into  two  groups;  H^,  Er,  Ey 
are  determined  through  scalar  function  of  B,  and  E  ,  Hr,  H0  are  determined  through 
scalar  function  of  A. 

Using  (2.89')  and  asymptotic  representation  P^[cos(Tr  -  9)1,  we  obtain  distant 
field  in  components  (in  CGS): 


dfrV  3iinS 


Here 
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(2.95) 


(2.96) 


For  calculation  of  field  in  layer  of  atmosphere  and  in  earth  it  is  necessary  to 
make  further  concretization  of  formulas  (2.95)  and  (2.96).  Let  us  introduce 
following  analytic  functions  of  two  arguments: 


IV. (**)  «, (*v) 

n';^)  I. 

U(*y)  «v(*y)  r 
’(,y)  !>,(*)  »',(*)  ’ 


[2.97) 
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where  Jv  and  nv  are  modified  Bessel  and  Neumann  functions  of  order  v,  while 

and  nj  are  their  derivatives  with  respect  to  argument.  Then  in  layer  of  atmosphere 

(V  rl); 

D,(«.r)|, 

Zj,r\ -  ,j  |  DJa 


(2.98) 


where 


%  ar  and  £  r-  MMi 


(2.99) 


A  (Ml 


*■  A(M) 


are  impedances  of  earth  to  spherical  waves  of  vertical  (?y)  and  horizontal  (*2Z) 
polarizations.  When  oQ  — ■  oo,  wave  number  of  layer  of  earth  kQ  -►  qd  and  —  0, 
while  1Z  -*  co. 

Coefficient  of  polarization  Hj(r)  of  normal  wave  of  number  J  determines  ratio 
of  intensities  of  components  (2.96)  to  components  (2.95).  It  is  calculated  during 
the  finding  of  eigenvalues  and  functions  of  operator  of  normal  waves. 

Taking  into  account  (2.98),  designating 

A{<v)  -  D,  [a,r)  -  Ov  («» , 

2, 


(2.100) 


0,  ( a',r )-  («,r) , 

E&7)-D,  J-D,  (•'/•) . 

-  n,  Z,DJaf) 

and  switching  to  system  of  MKS  units,  we  record  component  of  electromagnetic  field 
in  layer  of  atmosphere  in  the  following  form: 

W,  m IAJ \f  •~t*t J]  "J  {<**’>)  •  Dv*V)  . 

*  M 


(2.101) 
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'  /-« 

where  E  is  in  uv/m  and  H  is  in  pa/m, 

W  —  emissive  power,  in  kilowatts, 
x  —  wavelength  in  meters, 

A%  =  0.1829. 10”"  ‘  M',=  -4s-, 

Mj  »liy  ^  —  coefficients  of  excitation  of  wave  of  number  J, 

—  height  factor  at  point  of  radiation 


(2.102) 


(2.10J) 


and D*j(a',r')—  height  factors  at  point  of  recpetion, 

—  polarization  factor. 

Formulas  (2. 101) -(2. 102)  for  surface  of  earth  (r  =  a)  with  Hertz  doublet  also 
located  on  surface  of  earth  (b  -  a)  take  form: 
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(2. lot) 
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For  calculation  of  electromagnetic  field  under  surface  of  earth  or  ocean  we 


determine  Yj  and  Zj  in  interval  (0,  a): 


r) 


For  approximation  we  consider 


Atfcfl  ...  '*•" 


where 


H  =  a  -  r  is  distance  from  point  of  reception  to  surface  of  ocean  or  earth 


(2.106) 


8  is  attenuation  factor  of  flat  electromagnetic  waves  of  frequency  u)  penetrating 
medium  with  conductivity  aQ. 


Then  electromagnetic  field  at  depth  H  is  equal  to: 


W"  -JmJ  "  VA'  t) 
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(2.107) 


(2.108) 
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From  expressions  ( 2. 106)  (2. 108)  it  r'ollows  that  because  of  great  ground 
conductivity,  all  components  of  field  rapidly  diminish  with  submersion  of  point  of 
reception  in  water  per  the  same  exponential  law  with  attenuation  factor  6. 

§  7.  Generalization  of  Spectral  Theory  of  VLW  Propagation 
in  fcase  of^  bependence  of  Parameters  of  Medium 
on  Geographic  Coordinates  (Model  C) 

T^ethoa  of  Modulated  Normal  Waves) 

Solution  of  problem  of  VLW  propagation  on  the  basis  of  model  C  is  intimately 
connected  with  solution  of  problem  on  the  basis  of  model  B,  considered  in  §§  2-6 
of  this  chapter. 

According  to  §  1,  we  consider  that  model  C  is  characterized  by  layers  with 
radii  k  =  0,  1,  2,  . . .  N  -  1  weakly  depend  on  geographic  coordinates  0  and  qp. 

Let  us  record  this  fact  in  the  following  way: 

r*-rk  (pi,  w),  (2.109) 

where  p  is  small  parameter.  Coordinates  0  and  ep  enter  in  function  only  in 

combinations  p0  and  pep,  which  testifies  to  small  changes  of  r^  along  coordinates 

6  and  qp  during  the  period  of  wave. 

k 


We  also  consider  that  e. 


R  lr 

enn  and  depend  on  p0  and  pep,  that  is,  slowly 


"rrJ  c00 

change  from  coordinates  6  and  tp.  Under  these  conditions,  for  model  C  it  is  possible 

B 
A 


to  introduce  vector  function 


,  connected  with  components  of  field  according  to 


(2.25)-(2,J0)  and  determined  approximately  by  operational  equation 


(2.110) 


where  Lr  is  operator  generated  by  differential-matrix  expressions  (2.40)  and  also  by 
conditions  of  jumps  (2,42)-(2.43)  and  boundary  conditions  (2.44) -(2.45) .  Operator 


is  generated  by  differential  expression 


4. -tot  £(•"!"£•) 


1 


sin*#  <N* 


(2.111) 


and  conditions  of  limitedness  of  solutions  for  9  =  0  and  9  =  v. 

We  look  for  solution  of  equation  without  right  side  in  the  form 

m?). 


B  _ 

y  (r;  rt) 

A 

2  \r,  p?) 

(2.112) 


Putting  (2.112)-  in  (2.110),  for  I  =  0,  and  dividing  variables,  we  obtain 
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When  u  =  0  (model  3)  x  was  constant,  since  exact  splitting  of  equation  (2.113) 
in  two  took  place.  Now  with  small  |i  it  is  possible  to  consider  that  constant  is 
slowly  changing  function  of  coordinates  f?  and  d 

(2.U4) 

Then  we  obtain  from  (2.113)  equation  in  eigenvalues  of  x  of  operator  of  nomal 
waves  in  the  form 


L,  *  —  X( l*»  l*»)  * 


(2.113) 


Kigenvalues  of  Xj 

X*»  Xi»  /.*  ‘Aj  •••  (2.116) 

weakly  depend  on  coordinates  6  and  cp.  They  have  to  be  calculated  for  every  point 

k  k 

(0,  cp)  by  means  of  solution  of  problem  ( 2. 53) -(2.57)  >  in  which  parameters  e.£r,  tno, 
and  rk  are  determined  for  point  (0,  q?)  with  respect  to  assigned  functions 

•*  M  w);  »*?);  %  (A  M; 

To  every  eigenvalue  of  Xj  corresponds  eigenfunction 

Zi  ^1  y  =  0,  1,  2,  3  ...», 

\Zj  in  A  i*t)  l 

slowly  changing,  depending  upon  coordinate  9  and  cp. 

From  (2.113)  we  obtain  also  equation 


(2.117) 


determining,  function  of  phase  ^(<3,  4>) 
Free  normal  waves  (2,1 12)  have  form 


|4|-| 

Yj(r\^,  i»?) 

M/I  1 

• 

i/  n) 

If  I  l|*/(r;Aff)|  '  T/  (?118) 

They  were  called  [36,  37]  modulated  normal  waves.  We  borrowed  this  term  from 
work  of  .°>.  M.  Rytov  [42],  since  operator  of  modulated  normal  waves  it;  analogous  tc 
operator  of  modulated  normal  oscillations. 
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Applying  method  of  geometric  optics,  we  seek  solution  for  phase  factor  in  th 


form 


r,  -</,(Vri r 


where  Uj(0,  f)  are  slowly  changing  functions  of  0  and  q>. 

Then  for  functions  S,  and  U1  we  obtain  following  equations 


/- 0,14^ 


which  are  Fikonal  equations,  and 


expressing  law  of  conservation  of  energy. 

Thus  in  approximation  of  geometric  optics  modulated  normal  waves  have  form 

|*/l  JWM I  -M 

Ur  -1*^,1 1 


Equation  (2,120)  for  phase  Sj  of  normal  wave  determines  propagation  of  "radi 
beam"  in  two-dimensional  spherical  medium,  determined  by  coordinates  (0,  <p)  and 

characterized  by  local  "angular  refractive  indt 
v'-Xj  ,  which  changes  slowly  over  spherical  sun 


For  every  normal  wave  of  number  J  there  i 


particular  refractive  index 


Therefore 


radio  beams"  belonging  to  different  normal  wave 
pass  on  different  trajectories.  In  Fig.  5  are 
conditionally  shown  trajectories  of  "radio  beam, 
for  J  =  0,  1,  2,  ...,  connecting  point  of 
radiation  0(0  =  0)  with  point  of  reception 
P  (0,  qp).  Dotted  line  shows  geodesic  line. 


Quantity  S,  is  expressed  by  integral 


taken  on  radio  beam  connecting  point  of  radiation  0  with  point  of  reception  P*S 
determines  phase  of  normal  wave  of  number  j,  traveling  on  this  path.  Hero  u  :  I 
angular  element  of  arc,  taken  along  beam. 


U-i- ' 

u;«l) 

'  *fl 

- ■  -  ~  .*±*r**+  Vf  -  v.  ^ 

Distortion  of  "beam"  is  determined  by  usual  law  of  refraction,  which  follows 
from  (2.120).  Calculation  shov/^  that  for  first  types  of  normal  waves,  under  the 
condition  that  radio  route  does  not  pass  alonr  twilight  line,  distortion  of  beams 
is  so  insignificant  that  it  may  be  disregarded.  Therefore,  considering  that  "beam" 


spreads  along  geodesic  line,  we  determine  phase  Sj(j  =0,  1,  2)  by  integral 


=  jV-X,C, 


(M>)  A  , 


(2.12*0 


Puttirg  (2.124)  in  (2.121),  we  obtain  amplitude  factor  Uj  in  the  form: 


ij  at  /  -  X/  • 


Then  we  record  solution  of  (2.122)  in  the  form  of  modulated  normal  wave 

* 

\B>  ,***!'# 

U  “|z,  VmT‘  ■ 


! 


When  p.  =  0,  this  solution  becomes  asymptotic  representation  of  normal  waves  of 
model  B  already  obtained  in  §  B. 

Solution  (2.126)  is  applicable  both  to  "frozen"  and  also  to  si  wly  time-varying 
ionosphere.  In  last  case  Xj  may  be  represente  i  ir.  the  if  sum 

It  -fc  +  ty. 

where  Xj  is  time-averaged  value  of  Xj»  and  §Xj  is  variable  part. 

Accordingly,  in  formula  (2.126)  it  is  necessary  to  introduce 

Variation  of  parameters  of  ionosphere  can  lead  to  distortions  of  beams,  lowly 

varying  with  time.  We  shall  disregard  this  effect,  since  it  depends  on  member 
P  P 

d  v,/o<p  ,  which,  as  it  is  easy  to  show,  is  minute. 
iJ 

It  is  necessary  to  note  that  in  case  of  slowly  changing  ionosphere  equations 
(2.120)  and  (2,121)  are  applicable  only  under  condition  that  space  modulation  of 
ionosphere  has  shorter  periods  than  does  temporal  modulation,  that  is 

Xt  =  xM  w.i***)* 

where  p  is  small  parameter. 

In  our  case  periods  of  space  modulation  for  wavelength  X^  ~  20  Km  Is  lu-llj 
times  greater  than  Xq,  while  periods  of  temporal  modulation  are  greater  than  fieri  l 
T  =  ~  by  at  least  1000  times. 

-  |  J  _ _ 

•Linear  refractive  index  is  equal  to  - *  -  y  —If. 
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CHAPTER  III 


INTRODUCTION  OF  ADDITIONAL  INITIAL  DATA  OBTAINED  FROM  MEASUREMENTS 
OF  AMPLITUDE  AND  PHASE  OF  NEAR  FIELD  OF  VLW  AND  RESULTS 
OF  THEORY  ON  IONOSPHERE 

§  1.  Data  on  Near  Field 

O 

Since  parameters  of  medium  o>0  ^  and  r^  =  rQ  +  h^  (k  =  1,  2,  3,  ...  N  -  1)  are 
unknown,  then,  according  to  §  1,  Chapter  I,  we  solve  mixed  problem,  for  which  we 
introduce  additional  initial  data  on  near  field  Ej  and  Tj  at  points  (rQ,  ej)(J  = 

=  0,  1,  2,  ...  N),  located  on  surface  of  earth.  Such  data  were  obtained  in  great 
quantity  by  English  researchers  [31#  32,  33#  34,  35]#  who  measured  "instantaneous" 
distribution  of  amplitudes  of  near  field  along  earth's  surface  as  function  of 
distance  to  transmitter  [31]#  and  also  duirnal  variation  of  amplitude  and  phase  of 
components  of  field  at  fixed  points  at  distances  of  90  [32],  200  [33]#  and  535  km 
[34]  from  [VLW]  (CJXB )  station  with  call  sign  GBR  (16.0  kilocycles),  located  near  Rugby 
(52°  N,  1°  W).*  These  data  are  well-known,  since  investigation  of  near  field  has 
been  conducted  by  Cambridge  for  more  than  25  years  and,  starting  with  first  works  of 
Hollingworth  [52],  have  been  published  regularly. 

In  Figs.  4a  and  4b  are  given  typical  day  and  night  distributions  of  amplitude 
of  near  field  as  functions  of  distance. 

In  Table  2  are  given  distances,  where  there  are  maximum  and  minimum  field 
strengths  for  summer  period  (day). 


♦Signals  were  received  on  a  loop  in  two  positions:  in  vertical  plane  passing 
through  points  of  radiation  and  reception  (normal  component),  and  in  vertical  plane 
perpendicular  to  the  above-mentioned  (anomalous  component  of  field). 
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Table  2 


Max 

Min 

Max 

Min 

Max 

Min 

ad  (km) 

512 

155 

175 

224 

295 

490 

Dally  variation  of  amplitude  and  phase  of  normal  component  of  Held  for 
distance  of  90  km  is  given  in  Fig.  ?a  and  5b. 

Daily  variation  of  amplitude  and  phase  of  normal  component  of  field  for 
distance  of  555  km  is  given  in  Fig.  6a  and  6b. 

Analogous  dependencies  exist  for  anomalous  component  ( loop  perpendicular  to 
plane  of  propagation)  for  distances  of  90  km.  Comparison  of  normal  and  anomalous 
component  permits  us  to  separate  electromagnetic  field  reflected  from  ionosphere 
and  the  same  way  to  tie  i:  phase  in  Figs.  5  and  6  to  phase  of  direct  wave  [51-55]. 

It  is  necessary  to  note  the  exceptional  regularity  both  of  space  distribution. 
(Fig.  4)  so  also  of  daily  variation  of  amplitude  and  phase  of  near  field  (Figs.  5 
and  6),  which  are  repeated  from  day  to  day,  experiencing  slow  seasonal  changes. 
Exceptions  are  during  periods  of  flashes  of  radiation  on  sun  and  during  magnetic 
storms,  when  greater  or  lesser  considerable  deviations  from  above  mentioned  curves 
are  observed. 

§  2 .  Method  of  Solution  of  Reverse  Problem  of  Theory 
of  VLW  Propagation  for  Near  Field 

Data  on  near  field  depicted  in  Figs.  5  and  6  pertain  to  two  points  with  0  = 

=  0.014  radian  and  0  =  0.084  radian,  in  which  are  known  amplitudes  and  phases 

(E^ ,  tj)  and  (E0,  i'2)  at  any  time  of  day  and  year.  Thus  according  to  method  of 

solution  of  mixed  problem,  presented  in  §  1,  Chapter  I,  from  these  data  net  more 

than  4  parameters  of  ionosphere  can  be  determined.  Since  for  every  "k"  layer  of 

ionosphere  we  know  two  parameters  o)Q  ^  and  h^  =  r^  -  rQ,  then  data  of  Figs.  5  and  ( 

can  be  used  for  solution  of  problem  on  the  basis  of  two-layered  model  of  ionosphere, 

2  2 

for  which  parameters  h^,  and  cu^  h^  are  unknown.  Data  on  space  distribution 

of  near  field  (Fig.  4)  could  be  used  for  solution  of  problem  with  thinner  structure 
of  ionosphere,  but  they  are  incomplete.  Therefore  we  use  them  only  for  checking 
results  of  solution  based  on  two-layered  model  and  for  tying  in  phase. 
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Fig.  4. 
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According  to  diagram  of  Fig.  1,  we  have  to  solve  inverse  problem  of  theory  of 
waves  in  region  of  near  field,  which  is  reduced  to  solution  of  equations  (5)  and 
(6)  from  Chapter  I,  taking  following  form: 


|£)|  =  rood  F [P,  •;  (a,  0);  (fl,  •,);  a,  * •#,  (A,),  »,«(*,).  •{.» ,  «5. j,  A„A,) ; 


1>-  arg  F {P,  •,  (a,  0) ;  » «,oy) ; a,  7#;  mH,  *,„  ( A,),  (*»)»  «5. i»  •I.'*  *i>  *a) ;  ( 3  ■ • 1 

cont'd) 

J=U  2, 

where  P  -  moment 

a)  —  frequency 

(a,  0)  —  coordinates  of  Hertz  doublet 
(a,  0j)  —  coordinates  of  point  of  reception 
a  —  radius 

Fq  —  dielectric  constant  of  earth 

p  p 

“V  veff(hl)'  veff^h2^'  “o  1*  “o  2‘  hiJ  hp  -  data  on  first  and  second  layers  of 
ionosphere. 

Operation  F  designates  calculation  of  components  of  field  from  data  on  Hertz 

o 

doublet  and  medium.  Sought  quantities  for  system  of  equations  (5.1)  are  cuQ  ^  and 
2 

h^,  g ,  and  h2« 

Algorithm  of  calculation  of  components  of  field  in  the  form  of  spectrum  of 
normal  waves  represented  by  formulas  (2. 101) -( 2. 105) >  Chapter  II,  is  ineffective 
for  calculation  of  near  field  according  to  data  on  medium,  since  number  of  normal 
waves  which  have  to  be  considered  for  angular  distances  of  0  less  than  0.1  radian 
constitute  1.5-2  tens. 

Naturally  solution  of  inverse  problem  with  help  of  equations  (5.1) »  in  which 
operation  F  is  sum  of  15-20  members  depending  on  all  four  unknowns,  presents  great 

difficulties.  Therefore  we  switched 
from  spectral  representation  (2.101)  to 
beam,  in  which  every  component  of  field, 
for  instance  Er,  at  point  of  reception 
is  composed  of  series  of  beams 
experiencing  different  number  of 
reflections  from  first  and  second  layer 
of  ionosphere  and  from  earth. 

Thanks  to  anisotropy  of  ionospheric 
layers,  at  boundary  of  first  layer  and 

atmosphere  r  =  r^  =  a  +  h^  there  appear  two  refracted  beams,  ordinary  arid  extra¬ 
ordinary,  designated  by  letters  0  .  and  e,  .  On  second  boundary  of  ionosphere 

al  al 

r  =  r2  =  a  +  h2  incident  elliptically  polarized  beam  0d  or  ed  generates  two 


Fig.  7. 
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reflected  beams 


0  and  e  ,  remaining  in  first  layer,  and  two  refracted  beams  0 . 
rl  rl  a£ 


and  ea  ,  passing  into  second  layer  ( r0 ,  oo) 
a2  * 


Extraordinary  beams  experience  greater  than  ordinary;  therefore  basic  role  is 
played  by  ordinary  beams,  experiencing  small  number  of  reflections  and  refractions. 

In  Fig.  7  are  depicted  several  beams  remaining  ordinary  over  entire  section 
from  radiator  to  receiving  point. 

Direct  beam  for  emissive  power  of  1  kilowatt  has  intensity 

dq  is  len&th  of  path  of“  beam.  (3*2) 

Beam  that  has  experienced  one  reflection  from  lower  layer  of  ionosphere 

O  (5.2') 

Pi 


where  a  —  angle  of  incidence  on  earth 

P  —  angle  of  incidence  on  ionosphere 
A^  —  coefficient  of  focusing  of  ionosphere 
—  length  of  path  of  beam 

rBB  —  index  of  reflection  from  earth  of  vertically  polarized  beam 
R  —  index  of  reflection  from  ionosphere,  depending  on  parameters  of  ionosphere 

SB 

“0,1'  veff(hl)'  %  ^  rBB'  Rbb'  dSfined  in  work  [5]. 

Beam  experiencing  two  reflections  from  lower  layer  of  ionosphere, 

ep-t"  V<|[1  +  (3.2") 


where  RBr  and  RrB  are  conversion  factors  of  polarization  during  reflection  from 
ionosphere;  first  index  pertains  to  polarization  of  incident  ray,  the  second  to 
reflected; 


RBr  and  RpB  depend  on  angle  of  incidence  on  ionosphere  and  on  properties  of 
first  layer; 

rpr  and  rBB  are  indices  of  reflection  of  horizontally  and  vertically  polarized 

beams  from  earth;  remaining  designations  are  analogous  to  those  of  formula  (?.2). 

Beams  E^ ,  and  also  beams  E^),  E^ . . . ,  refracted  by  first  boundary 

of  ionosphere  and  reflected  from  second  are  determined  similarly. 

Besides  ordinary  E)  ' ,  determined  by  propagation  along  directions  0.  ,  0  in 
i  al  rl 

layer  (r^,  Tg),  extraordinary  beams,  propagating  in  directions  ed  ,  er  ,  and  also 


mixed  beams  0 .  ,  e..  ,  and  e . 

d1  r^  d3 


,  0  are  possible.  All  of  them  experience  one 

*  A 
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mam 
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reflection  from  layer  r2  in  passing  through  boundary  of  layer  r  =  r^.  Besides 

/  p) 

ordinary  beam-  EX  ' ,  determined  by  direction  0,  ,  0_  ,  0.  ,  0  ,  there  are  15 

d.  Q1  ri  Q1  ri 

additional  beams  of  mixed  type  and  purely  extraordinary  beams  e,  ,  e„  ,  e ,  ,  e  , 

al  rl  al  rl 

experiencing  two  reflections  from  ionosphere. 

Full  field  is  composed  of  infinite  number  of  beams  experiencing  one,  two, 
three,  etc.,  reflections  from  ionosphere 

£  =  J3 


»-o 


(3-3) 


All  of  them  may  be  calculated  from  laws  of  geometric  optics,  if  matrix  of 
indices  of  reflection  from  earth  is  known 


(3.4) 


and  matrix  of  indices  of  reflection  and  refraction  of  beams  on  boundaries  of  layers 
r  =  r^  and  r  =  r2. 

For  lower  boundary  of  ionosphere  r  =  r^  there  have  to  be  given  indices  of 
reflection 

i-M,  *V\ 


#« 


and  refractive  indices 


D«*»  D<" 

n  •  n 


(3.5) 


(3.6) 


Analogously,  for  second  boundary  of  ionosphere  r  =  r2  there  have  to  be  known 
matrices  of  indices  of  reflection 


(3.7) 


and  refraction 


M 


Rm 

Aw, 

* 

P*21 

*'r»» 

* 

Dm 

tf. 

n<*' 

LJt%  • 

1 t 

(3.8) 


For  calculation  of  near  field  of  VLW  with  accuracy  to  5#  it  is  sufficient  to 
consider  only  first  five  beams  depicted  in  Fig.  7: 

£*£,  +  £Y>  4- £?!  +  £?> +  £*>>.  (5.9) 

For  small  distances  (to  100  km)  beams  Eq  and  dominate.  For  large  distances 

(300-600  km)  beams  EQ  and  E ^  dominate. 

At  distances  over  1000  km  beam  representation  (3.3)  is  impermissible,  since  for 
beams  with  small  number  of  reflections  caustic  of  beams  near  surface  of  earth. 
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Therefore  for  sliding  beams  it  is  necessary  to  consider  phenomenon  of  diffraction. 

As  will  be  seen  from  Chapters  IV  and  V,  the  main  part  of  distant  field  of  VLW 
is  created  by  beams  not  touching  earth's  surface.  They  move  on  ricochetting 
trajectories  (adhesion  effect).  Thu-,  distant  field  of  VLW  on  earth's  surface  is  due 
only  to  seepage  of  waves  beneath  surface  of  caustic  of  beam  and  therefore  carries 
purely  diffractive  character. 

During  calculation  of  field  strength  of  beams  EQ,  E^1),  E ^ ...  it  was 
assumed  that  coefficients  (3.4)-(3.8)  can  be  approximated  by  corresponding 
coefficients  for  flat* medium,  but  for  calculation  of  length  of  path  of  beams  Dq, 
D^1),  ...  sphericity  of  layers  of  ionosphere  and  earth  was  considered.  Dis¬ 

regard  for  sphericity  evokes  impermissible  error  in  phase  of  beams,  since  with 
distance  of  500  km  flat  and  spherical  models  give  difference  in  the  course  of 
ionospheric  beams  E  of  2-3  km.  First  model  can  be  applied  for  calculation  of  fields 
only  to  distance  of  not  more  than  300  km. 


2  2 

Determination  of  quantities  1,  a>0  2,  h^,  and  h2  from  equations  (3.1)  was 
carried  out  with  help  of  the  presentation  of  field  in  the  iorm  (3.9)  by  graphoanalytic 
method,  which  will  bo  described  in  another  place.  Here  we  shall  give  only  results 
of  calculation  for  vernal  equinox,  which  with  small  changes  can  be  applied  to  whole 
spring  and  summer  seasons  and  partially  to  autumn. 


2  2 

On  Fig.  8a  and  8b  are  depicted  distributions  u>Q  ^  and  u>0  2  obtained  from 
calculation,  depending  upon  height  of  layer  h  at  noon  (8a)  and  at  night  (8b). 
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Corresponding  values  of  v&pf  were  determined  by  Nicolet  graph,  (Chapter  II, 

§  1)  and  plotted  in  Fig,  8  on  corresponding  steps  of  curve  cu^h) . 

Check  of  two-layered  model  of  ionosphere  for  daytime  distribution  of  field 
(Fig.  4)  shows  what  6  points  of  maxima  and  minima  of  interference  curve  (3.9)*  marked 

in  Table  II,  conform  well  to  diagram 
depicted  in  Fig.  8a.  This  indicates 
that  what  actually  existing  distribution 
of  ionization  in  daytime  is  character¬ 
ized  by  two  regions  of  height  h^  and  h2 
with  sharp  gradients  of  ionization. 
Englishmen  [48]  call  these  layers 
Dp(h  =  h1)  and  DQ(h  =  h2) . 

For  nighttime  (Fig.  8b)  we  do  not 
have  reliable  graphs  of  dependence  of 
field  E  on  distance  similar  to  Fig.  4b, 
but  it  is  possible  to  assume  that  two-stage  picture  of  ionization  depicted  in  Fig.  8b 

p 

is  result  of  idealization  of  problem  and  that  in  reality  u>Q(h)  is  characterized  by 
smooth  drop,  the  "trail"  of  ionization  of  layer  E  at  night.  It  is  interesting  to 

p 

note  that  jumps  of  <aQ  on  boundary  h^  by  day  and  night  are  approximately  identical. 

Data  on  daily  variation  of  phase  (Figs.  5  and  6)  permit  us  to  find  daily 
variation  of  heights  h^  and  h^,  depicted  in  Fig.  9.  As  can  be  seen  from  Fig.  9, 
lower  bound  of  layer  of  ionosphere  h^  follows  geometric  shadow  of  earth  up  to  night 
layer  E^.  Second  layer  h2  is  connected  with  height  of  sun  \c  by  empirically  found 
dependence 

h  A,  -j-  A  In  see  Xff  (5.10) 

where  h0  «  70  km  and  A  is  slowly  changed  during  the  course  of  the  year  from  3  to  7 
[31]. 

Formula  (3.10)  is  valid  to  xc  =  ®5°.  For  distant  field,  as  will  be  seen  later, 
basic  role  is  played  by  lower  layer  h  =  h1. 
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CHAPTER  IV 


BASIC  CHARACTERISTICS  OF  NORMAL  WAVES  OF  FIRST  NUMBERS 
§  1.  Criteria  of  Similarity 

Having  completed  calculations  made  in  Chapter  III,  we  know  all  parameters  of 
medium.  Therefore  it  is  possible  to  turn  to  last  operation  of  solution  of  mixed 
problem  per  diagram  of  Fig.  1,  calculation  of  distant  fVLW]  (CflB)  field  by  spectral 
formulas  ( 2. 101) -(2. 108)  (model  B)  and  (2.127),  (2.128)  (model  C).  In  these  spectral 
formulas  enter  following  characteristics  of  normal  waves: 


1)  wave  numbers  v.  =  a.  +  ip., 

J  J  J  /T  2 

2)  coefficients  of  excitation  nj  =  XV/  ^-Vy  i 


3)  polarization  factors 


ZAr)  >r)~~Z,Dv  (o,r) 

*i"TT7T - 1 - 1  - - 

ves  Y,( r) . 


(*.l) 


(4.2) 


4)  forms  of  normal  waves  Yj(r). 

In  model  B  tiny  are  constant,  and  in  model  C  the  change  slowly,  depending  upon 
geographic  coordinates  9  and  cp.  Calculation  of  Vj,  rij,  >tj  and  Y  for  both  models  is 
reduced  to  finding  eigenvalues  of  Xj  and  eigenvectors  |Yj|Zj|  of  operator  of  normal 
waves,  described  by  expressions  (2.53)-(2.57)  for  particular  case  of  four-layer 
medium  k  =  0,  1,  2,  3,  where  layer  k  =  0  is  earth,  k  =  1  is  atmosphere,  and  k  =  2 
and  k  «  3  are  layers  of  ionosphere. 

Thus,  Vj,  nj,  Hj,  Yj  are  functions  of  the  following  9  dimensionless  parameters, 
characterizing  medium  at  given  point  9,  <p: 

1)  radius  of  earth  k^a  «  kj_ro»  where  k^  =  p 

.  i  m 

2)  dielectric  constant  of  earth  eQ  =  e0  +  ie0. 
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3)  height  lower  boundary  of  first  xuyer  k-jh^  =  ki(r1  -  rQ) , 

4)  square  of  critical  frequency  of  first  layer  u>q  ^/a>‘  , 

5)  number  of  collisions  in  first  layer 

6)  height  of  lower  boundary  of  second  layer  =  k^Tg  "  tq)  > 

7)  square  of  critical  frequency  of  second  layer  g/cb , 

fi)  number  of  collisions  in  second  layer 

9)  gyromagnetic  frequency  o^/ui. 

From  expressions  for  operator  of  normal  waves  (2.53)-(2.57)>  and  also  (4.1) 
and  (4.2),  it  follows  thas  Vy  riy  Wj  and  Yj  are  analytic  functions  of  these 
parameters.  Contemporary  mathematical  means  do  not  permit  us  to  find  these  functions 
in  analytic  form;  therefore  we  used  ordinary,  well  approbated,  direct  mathematical 
methods  of  finding  eigenvalues  and  functions,  which  give  numerical  values  of  Kj  and 
|  Yj( r)/Zj(r) |  for  given  numerical  values  of  parameters  of  medium  k^a,  e0J  k^h.^;  k^hg, 
etc.  Eigenvalues  of  Xj»  and  this  means  of  Vj  also,  were  calculated  with  accuracy 
to  5  significant  digits,  while  Yj(r)  and  Zj(r)  were  calculated  to  3  significant 
digits.  Coefficients  of  excitation  and  polarization  factors  were  found  from 
calculated  v j  and  |Yj/Zj|  on  tne  basis  of  formula  (4.1)  and  (4.2)  with  accuracy  of 
3  significant  digits. 

Exact  calculation  of  Vj,  nj,  Kj  and  Y^  over  wide  range  of  changes  of  nine  shown 
parameters  of  medium  by  direct  algebraic  methods  involves  a  great  deal  of  work, 
which  scarcely  need  be  done  for  clarification  of  basic  regularities  of  distant  VLW 
field.  For  reduction  of  volume  of  calculation  operations  we  adhered  to  the  following, 
which  were  checked  by  calculations: 

1.  At  distances  of  more  than  2000  km  basic  contribution  to  spectral  formulas 
of  models  B  and  C  in  range  of  frequencies  of  10-20  kilocycles  is  made  by  first  three 
normal  waves,  called  quasi-TH^,  quasi-THg,  and  quasi-TE^.  Remaining  types  of  normal 
waves  may  be  omitted  from  calculation. 

2.  For  shown  types  of  normal  waves  basic  role  in  waveguide  process  is  played 

by  lower  bound  of  first  layer  of  ionosphere  h^  =  -  rQ  and  surface  of  earth, 

whereas  for  high  types  of  waves  the  role  of  waveguide  is  played  by  boundary  of  second 
layer  of  ionosphere  and  surface  of  earth  (see  Fig.  11).  Therefore  characteristic 
criteria  of  similarity  for  first  types  of  waves  will  be: 

•ai  TUi 

M ;  «,  =*  •',+*  «V  *,*, ;  Ag- ;  ;  — . 
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5.  For  frequencies  in  10-20  kilocycle  range  number  of  characteristic  criteria 
can  be  further  reduced,  since 

.  .  4w# 

a)  in  virtue  of  high  conductivity  of  soil  and  sea  water  — 

• 

i 

and  e  q  is  immaterial 

b)  under  the  condition  that  frequency  o>(~10^)  is  almost  one  order  of  magnitude 

6  “r 

less  than  v  ff(~10  ),  components  of  tensor  ||e||  depend  only  on  two  parameters  and 


vr  ■ 


1+i 


•  cost, 


sin  t , 


where 


m 


cost, 


•ft 


V/ 


Thus  in  working  range  of  frequencies  m  wave  numbers  Vj  are  functions  of  five 
criteria: 


here  <ur  ^/u)  and  pertain  to  lower  layer  of  ionosphere.  On  these  criteria  of 
similarity  of  medium  also  depends  remaining  characteristics  of  normal  waves  nj,  Hj» 
and  Yj. 

For  study  of  field  of  VLW  at  the  surface  of  earth,  where  Yj  ■  1,  it  is  necessary 
to  know  only  three  characteristics  of  normal  waves  Vj,  nj,  and  y.y  as  functions  of 
5  above-indicated  criteria  of  similarity  of  medium  in  ranges  corresponding  to 
possible  values  of  parameters  of  medium  on  route  at  different  times  of  day  and  year. 

Before  describing  results  obtained,  let  us  divide  discrete  spectrum  of  normal 
waves  into  two  series  of  normal  waves,  called  quasi-THj  (j  =  0,  1,  ?,  7,  ...  arid 
quasi-TEj  (j  *-  1,  2,  5,  ...).  By  definition,  first  types  of  waves,  for  t  -*  0, 
continuously  pass  into  ordinary  waves  of  type  THj  (J  =  0,  1,  2,  ...)  with  components 
Hjp,  Er,  Eq,  while  second  become  ordinary  TEj  waves  with  components  E^,  Hr  and  H@. 

THj  and  TEj  waves  determine  electromagnetic  field  in  Watson  model,  where 
magnetic  field  of  earth  is  absent  (t  =  0).  Here  vertical  Hertz  doublet  excites  only 
TH -waves . 
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In  case  of  t  /  0  waves  of  quasi-THj  type  with  wave  numbers  (J  =  C,  1,  2,  3, 
. ..),  besides  basic  components  H^,  Er  and  Eg,  contain  components  E^,  Hr,.  and  H0. 
Ratio  of  two  groups  of  components  is  determined  by  magnitude  of  polarization  factor 
hJH*  Analogously,  quasi-TEj  waves  with  wave  numbers  VjE,  besides  basic  components 
E^,  Hr#  and  H0,  contain  components  H^,  Er,  and  E0  in  ratio  determined  by  quantity 
1/kjE.  With  sufficiently  large  t,  HjE  and  become  comparable  to  1;  therefore 
criterion  of  domination  of  components  E^,  Hr,  and  H0  in  quasl-TEj  waves  and  of 
components  H^,  Er,  E0  in  quasi-THj  waves  cannot  be  used  for  discerning  THj  and  TEj 
waves.  For  their  identification  during  calculations  we  continuously  decreased  t  and 
0. 


It  is  natural  that  when  x  /  0,  vertical  Hertz  doublet  excites  both  qvmsl-THj 
and  quasi-TEj  waves.  Coefficients  of  excitation  of  these  waves  subsequently  are 
designated  njH  and  HjE  respectively. 

§  2.  Dependence  of  Charactex'istlcs  of  Normal 
Waves  on  Wave  Number 


o 


In  Fig.  10  are  depicted  spectra  of  normal  waves  of  type  THj  for  j  =  0,  1,  2,  3, 
4,  5  and  6  for  following  parameters  of  medium:  a  =  6370  km,  crQ  =  oo,  h^  =  70  km, 

0)  A  U)  0  p. 

hg  -  75  km,  -  0.15j  ■—*—  =  1;  tg  t1  *=  0.4;  tg  Tg  =  0.8,  and  u>  =  lcr,  which 

corresponds  to  daytime  conditions  of  propagation  above  marine  surface.  On  abscissa 
are  plotted  ct j  =  Re  (vj),  while  ordinate  are  plotted,  upwards,  mod  n^  and,  downwards, 
3 j  »  Im  ( v j ) . 

In  Fig.  11  are  depicted  distributions  of  amplitude  and  phase  q>  of  component 
ErJ  respect  to  height  for  J  =  1,  2,  and  4  (TH1,  THg,  and  TH^). 

On  the  basis  of  Figs.  10  and  11  the  following  conclusions  can  be  made. 

1)  Real  parts  of  wave  numbers  ctj  =  Re  (tj)  monotonically  decrease  with  growth 
of  wave  number.  The  biggest  value  of  a0  for  TH0  wave  is  always  smaller  than 
magnitude  of  k^r^,  which  is  due  to  positive  value  of  angle  of  slip  y  with  incidence 
of  normal  waves  on  layer  of  ionosphere: 

COIty*  fifty)/*, r,.  (4.3) 

2)  Period  of  modulation  of  amplitude  of  normal  wave  (Fig.  11)  with  respect  to 
front  decreases  with  increase  of  wave  number  j. 

3)  Coefficients  of  excitation  nj  of  normal  waves  sharply  decrease  for  first 
types  of  waves,  for  which  2  k^a,  which  is  caused  by  low  field  strength  of  these 
waves  near  earth's  surface  (adhesion  effect). 


-62- 


4)  Attenuation  factors  Pj  =  Io  (yj)  of  normal  waves  monotonically  increase  with 
increase  of  number  j.  Exception  is  THq  wave,  which  experiences  large  degree  of 
damping.  THq  wave  strikes  ionosphere  at  an  angle  close  to  quasi-Brewster  angle; 
therefore  it  is  reflected  weakly  from  ionosphere. 

5)  From  Fig.  11  one  mey  see  that  first  types  of  waves  are  reflected  basically 
from  lower  boundary  of  first  layer,  and  high  types  of  waves,  starting  with  TH^,  are 
reflected  from  second  layer. 
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Fig.  10. 


Numerical  values  of  nj  are  multiplied  by  2/?*101®, 

Of  special  interest  is  THq  wave,  which  is  analog  of  cable  wave  in  two-conductor 
waveguide  systems.  In  flat  models  of  media  it  is  considered  the  basic  wave  [27,  28, 
29].  During  calculation  of  sphericity  of  earth,  cable  wave  loses  dominating  value, 
because  its  field  is  concentrated  near  concave  surface  of  ionosphere  and  is  weakened 
greatly  near  earth  (adhesion  effect);  coefficient  of  excitation  nQ  <  <  nj  (.1  /  0). 
Besides,  in  daytime  cable  wave  possesses  experiences  greater  degree  of  damping  from 
incidence  at  angle  close  to  Brewster  angle. 

From  Figs.  10  and  11  it  follows  that  in  daytime  at  distances  over  1000  km  it 
is  necessary  to  consider  only  quasi-TH^  and  quasi-THg  waves. 


For  nighttime,  dependencies  of  Vy  n^ ,  Mj  and  Yj  on  number  J  have  analogous 
form.  Number  of  propagating  normal  waves  of  quasi-TH  type  is  increased  to  8 ; 
therefore  Aa^  =  aj  -  decrease.  Besides,  all  atj  are  displaced  to  the  right; 

therefore  larger  number  of  first  types  of  normal  waves  is  subjected  to  adhesion 
effect,  showing  up  as  decrease  of  coefficients  of  excitation  nj.  At  night  qausi-TEj 
waves  start  to  play  a  role.  Their  wave  numbers  ajE  =  Re(vjE)(J  =  1 t  2,  3,  ...)  also 
monotonically  decrease  with  growth  of  number  J,  starting  from  first  TE^  wave,  and 
attenuation  factors  monotonically  increase. 

Increase  of  excitation  frequency  leads  to  change  of  spectra  (Fig.  10),  to 
analogous  increase  of  height  of  layer  h^.  Number  of  THj  and  TEj  waves  is  increased, 
that  is,  AOj  decrease,  and  all  are  displaced  to  the  right.  Therefore  with 
increase  of  u>,  an  ever  greater  number  of  first  types  of  normal  waves  falls  in 
interval  k^a  <  a  <  k^(a  +  h1),  where  adhesion  effect  appears.  Thus  they  cease 
playing  role  in  spectrum  of  distant  field,  and  their  place  is  taken  by  higher  types 
of  waves. 

Calculation  shows  that  in  range  of  frequencies  from  10  to  20  kilocycles  and 
with  change  of  height  of  layer  h1  from  70  to  90  km  at  distances  of  more  than  2000  km, 
as  a  rule,  basic  role  in  distant  field  is  played  at  night  by  three  waves:  quasi«TH1, 
quasi-TH2,  and  quasi-TE^.1  Exception  to  this  occurs  in  certain  sections  of  route, 
where,  as  a  result  of  interference,  there  occurs  mutual  extinguishing  or  quasi-TH^ 
and  quasi-TH2  waves  (0  =  eminmin'  Chapter  V,  §  2).  At  these  singular  points  it  is 
necessary  to  consider  higher  types  of  waves. 

Subsequently  we  shall  limit  consideration  to  but  three  types  of  normal  waves, 

TH^,  TH2,  and  TE^,  since  above-indicated  sections  of  route  are  unfit  for  communication 
and  navigation. 

§  3.  Dependence  of  Characteristics  of  Normal 
Waves  on  Properties  of  Ionosphere 

Properties  of  ionosphere  are  determined  by  three  criteria  of  similarity  kjh^, 
u>r  ^/uj  and  In  Figs.  12,  13,  and  14  is  given  complex  plane  v  =  a  +  iP,  on  which 

are  plotted  curves  of  v^,  v^,  and  as  functions  of  parameter  o>  ^/b  of  first 
layer  of  ionosphere  for  two  different  values  of  parameters  t  and  k^h,  with  u>  =  10  : 

a)  tgtj-0,4  (vjy-2.101), kxhx  -  22(A,=70»m),  A.-75-,  *«*}=  10*. 
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Igt,  -  0,4  (i«|>  -  2.10’), Mi  -  28  (A,  -88,5m).  h,  -  93- ,  *<J » - 1  O’ , 


tgt  -  10K''=  8.10s),  Mi  =  22(4,-70.-),  At-76«,,^)-  5.10s, 


d) 


-^  =  1,5; 

tgt  -  10  -  8.10s),  M,*=28 (A, =88,5.-),  A,— 93  — , —  5.10s, 

•*»/•- 1.5* 


Remaining  parameters  of  medium  had  values  k^a  =  2001.2;  aQ  =  oo. 

Increase  of  cd  ./od  along  every  curve  (Figs.  12-14)  occurs  because  of  increase 
r,  i 


of  electron  concentration  N  '  ' 

e 


of  first  layer  of  ionosphere,  that  is,  of  quantity 

<J0 


OJ, 


■0  1  (In  Figs.  12-14  are  plotted  values  of  — . 

Real  parts  of  v^H,  v2H,  and  monotonically  grow  with  decrease  of  ciy^co. 

Changes  of  dampings  P^,  ^2H’  and  ^1E  as  factions  40 r are  more  complicated. 
The  most  simple  is  dependence  of  P1£  on  co^/a  for  quasi-TE^  wave  (Fig.  14).  With 
decrease  of  —  from  oo  (ideally  conducting  ionosphere),  damping  P1E  =  Im  (v.jV) 
monotonically  increases  from  0.  Here  real  part  of  wave  number  a..  =  Re  (r.,,)  is 

increased  proportionately,  so  that  line  becomes  straight.  Increase  of  P1E 

is  explained  by  weakening  of  reflective  properties  of  h^  layer  of  ionosphere. 

For  TH^  and  TH^  waves  monotonic  build-up  of  PjE  with  decrease  of  u>r/cD  is 
disturbed  in  range  of  values  of  cd^/co  approximately  from  10  to  0.5,  where  there  is 
observed  a  quite  sharp  fall  P1E  and  P2H  with  decrease  of  0^/10. 

When  u^co  =  10  to  20,  maximum  damping  occurs.  Appearance  of  maximum  damping 
is  due  to  fact  that  normal  waves,  quasi-TH^  and  quasi-THg,  with  predominance  of 
vertical  polarization  (n^  and  are  small) ,  strike  layer  of  ionosphere  r  =  at 
angles  close  to  quasi-Brewster  for  values  of  a>r/u>  =  10  to  20.  Therefore  in  this 
range  of  values  of  cd^/oo  reflective  qualities  of  ionosphere  are  weakened.  Subse¬ 
quently,  for  brevity,  we  shall  call  this  phenomenon  the  Brewster  effect. 


CD 


r,  1. 


When  ^ ti s  are  considerably  smaller  than  quasi-Brewster  values,  there  is 


co 


restored  normal  increase  of  P,„  with  decrease  of 

Jn 


nl 

(0 


,  caused  by  weakening  of 


reflective  qualities  of  ionosphere.  Between  this  normal  region  and  region  of 

co  1 

quasi-Brewster  maximum,  in  interval  5  >  — >  0.?,  damping  is  minimum,  wherein  lie 
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values  of  v,u  and  v„„  corresponding  to  normal  conditions  of  lower  layers  of 

1H  cn 

ionosphere. 


Fig.  12. 


I 

i 

| 


For  values  of  o>  ../id  less  than  0.2  lower  layer  of  ionosphere  is  so  greatly 
weakened  that  basic  role  is  played  by  higher  disposed  second  layer  of  ionosphere. 

From  comparison  of  curves  (Figs.  12-14),  with  tg  t  =  0.4  and  tg  t  ■  10,  it 
follows  that  magnetic  field  of  earth  starts  to  show  up  with  large  values  of  t,  when 
v  ff  is  less  than  u^.  Here  attenuation  factors  decrease.  Influence  of  magnetic 
field  of  earth  is  greatest  for  wave  quas.i-TH1  waves,  less  for  quasi-TH2  waves,  and 
practically  negligible  for  quasi-TE^  waves.  Decrease  of  and  due  to  magnetic 
field  of  earth  is  observed  over  wide  range  of  value  of  u>r;j/cD  <  5  and  is  explained  by 
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A. 


weakening  of  effect  of  Brewster  angle  as  a  result  of  appearance  of  magnetic 
anisotropy  of  layer.  In  region  of  large  values  of  a>rl/cu  >  5,  where  Brewster  effect 
is  maximum,  influence  of  magnetic  anisotropy  is  small. 


9<i 


From  Figs.  12-14  can  be  made  also  conclusions  about  altitude  effect  of 


cu 


ionosphere  on  Re  (vj)and  Jm  (v^).  When  —  >  0.2-0. 3,  real  parts  of  wave  numbers 
a^,  agH,  and  a1E  grow  with  h^,  just  as  for  all  other  Oj  (Fig.  10).  Here  speed  of 
3alH 


displacement  is  least;  therefore  intervals  between  and 


and  also  between 
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a^H  and  a^,  decrease  with  growth  of  h.  This  regularity,  shown  earlier  during 
discussion  of  spectra  (Fig.  10),  very  significantly  changes  picture  of  distant  field 
during  transition  from  day  to  night.  Thickening  of  spectrum  with  increase  of  width 
of  waveguide  h^,  or  of  excitation  frequency  (criterion  k^h^),  is  general  dependence 
for  all  hollow  waveguide  systems. 

Attenuation  factors  P^,  ?2H'  ^1E  ^rop  w^h  increase  of  h. 


§  4.  Dependence  of  Characteristics  of  Normal  Waves 
on  Properties'  of  EarfcK 


47rj0 

Properties  of  earth  are  determined  by  parameters  k^a  and  — — .  Influence 
earth's  curvature  on  wave  numbers  in  range  rQ  -  6 350-6380  km  is  well  within 
following  linear  dependence: 


of 


«;('«)  _  •>(*) 
r%  a 


(4.4) 


where  a  -  6370  km. 

For  high  types  of  waves  this  dependence  is  satisfied  in  wider  limits  than  for 
first  types  of  waves,  where  this  dependence  is  distorted  by  effect  of  "adhesion." 

With  increase  of  radius  of  earth  rQ,  effect  of  "adhesion"  is  less,  and  formula  (4.4) 
is  satisfied  with  large  degree  of  accuracy.  Conversion  of  from  plane  case 
(a  -*oo)  to  spherical  with  a  =  6370  km  by  formula  (4.4)  for  first  types  of  waves  is 
impermissible,  which  indicates  unfitness  of  flat  models  for  calculation  of  distant 
field  of  VLW.  When  a  =  6370  km,  form  of  TH1  wave,  determining  distant  field,  differs 
considerably  from  TH^  wave  of  flat  case  (a  =00). 

In  range  of  radii  of  earth  of  6350-6380  km  it  is  possible  to  consider  that 
forms  of  waves,  coefficients  of  excitation  n^,  damping  Pj,  and  polarization  Hj 
remain  constant  under  the  condition  that  other  parameters  of  medium  k^r^,  t  etc., 
are  not  changed. 


More  essential  influence  on  characteristics  of  normal  waves  is  rendered  by 
conductivity  of  "soil."  In  Fig.  15  is  depicted  dependence  of  attenuation  factors 
of  quasi-TH^  and  quasi-THg  waves  on  resistance  1 /oq  for  day  (h^  =  70  km,  = 

=  2.  10  ,  etc.)  and  night  (h  =  88.5  km,  =  8.10  ,  etc}  cases  when  a)  =  10J.  As  can 

be  seen  from  graph,  in  range  of  "soil"  resistances  from  sea  water  to  very  dry  soil 
attenuation  factor  of  TH^  wave  is  increased  by  0.7,  while  that  of  TH^  wave  is 
Increased  by  0. 9-1.0. 
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Pig.  15. 

For  approximate  calculation  of  influence  of  conductivity  for  all  types  of  waves 
it  is  possible  to  use  formulas: 

PA)  •  M°°)  +  A?A). 

*/*•)  *=  */*)+  &i/Vi 

where  Vj(oo)  =  aj(a°)  +  iPj(ac)  is  wave  number  of  normal  wave  J  for  oQ  -  oo;  quantities 
Aaj(°o)  and  APj(°0)  are  Presentecl  ln  Pig.  15;  for  TH.^  wave  is  plotted  separately, 

since  Influence  of  Oq  on  TH^  wave  is  somewhat  less  than  for  high  types  of  waves, 
owing  to  effect  of  adhesion. 

X2V^10“ 


Fig.  16. 
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On  the  same  Fig.  15  is  given  dependence  of  increase  of  Acij  with  impairment  of 
ground  conductivity.  Aa1H  sma-'-1  anJ  Is  not  shown  in  Fig.  15. 

In  Fig.  16  is  given  dependence  of  coefficients  of  excitation  n^  for  quasi-TH^ 
and  quasi-THg  waves  on  resistivity  of  soil  l/a0.  As  can  be  seen  from  graphs  of 
Fig.  16 ,  influence  of  1/cJq  on  nj  is  insignificant.  Just  as  insignificant  is  influence 
of  resistance  of  soil  on  coefficients  of  polarisation  k^. 

§  5.  Dependence  of  Characteristics  of  Normal  Waves 
on  Parameters  of  Radiator 


Radiating  device  in  the  form  of  vertical  Hertz  doublet  is  characterized  by 
emissive  power  W,  coordinates  of  dipole,  and  frequency  u>. 

Emissive  power,  as  follows  from  formulas  (2.101)  and  (2.103),  determines 
amplitudes  of  normal  waves,  which  are  proportional  to  s'quare  root  of  emissive  power. 
Characteristics  of  normal  waves  Vy  v.y  nj,  and  Yj  do  not  depend  on  emissive  power. 

Position  of  dipole,  determined  by  height  of  dipole  above  earth  b  -  r  ,  does  not 
affect  Vj  and  Hj,  but  can  strongly  change  coefficients  of  excitation  nj(0)  at  point 


Frequency  o>  of  radiator  essentially  affects  all  characteristics  of  normal 
waves,  since  it  enters  in  characteristic  criteria  k^a,  UwOq/x,  oo^/cd,  and  kjh^. 
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General  character  of  regularities  of  change  Oj  and  Bj  for  TH^,  TH2,  and  Tl^ 
waves  from  frequency  o>  may  be  established  from  Figs.  12-14.  With  increase  of  u), 
criterion  k.^  is  increased,  while  tu decreases,  where,  according  to  Figs.  12-14, 
all  otj  increase  and  B 
maximum  and  increase  for  region  of  small  oj^/cu.  Since  influence  of  frequency  to  on 
characteristics  of  normal  waves  is  complicated  effect,  in  which  participate  4 
parameters  of  medium,  then  for  convenience  of  calculation  in  Figs.  17-22  are  given 
to  dependencies  of  the  following  quantities  on  frequency: 

-  *1*!;?;;  mod n/,  arg*/  mod*/  arg*y. 

On  every  figure  are  given  two  curves  for  "day"  and  "night"  conditions.  They 
are  basic  working  graphs  for  calculations  of  amplitude  and  phase  of  distant  field 

(Chapter  V).  For  "day"  curves  are  taken  parameters  of  medium 

o  =  6370  KM,  A,  =  69,73  km:  «>•».,  =  0,3-10**,  »<)>=2-10T, 

A,  =  75  km,  •|.j  =  2-10,,!  v<J >  =  107,  ^  =  8-10*. 

For  "night"  curves  are  taken  parameters  of  medium 


1H  and  B2H  decrease,  if  to^/co  lies  in  tail  of  quasi-Brewster 


0  =  6370  km,  A,  =86.5  km,  ®J.i 
*<]»  =  8- 10s,  A,  =90  km,  <2  =  1,5- 10”,  vP)  « 


=  0,3- 10”; 


5-10\  •U/y^S-lO*. 
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In  addition,  in  Fig.  17  ’s  given  dependence  =  a^H  (86.5)  -  a  j,  (70). 

Magnitudes  of  a 2H  and  a1E  are  calculated  by  the  formulas:  o^h  (70)  =  (70)  - 

"  (alH  '  a0H^  a2H  (S6*5)  "  aiH  (7°)  +  AalH  '  (alH  '  a2H^  etc»  Where  alH  (7°)  = 

=  1611.4;  2004.0;  2604.4,  respectively  for  oj  =  0.805;  1.000;  1.290  x  0.942*105. 

Parameters  of  waves  for  intermediate  heights  are  interpolated  per  linear  law. 

For  calculation  of  conductivity  of  earth,  Figs.  17-19  for  cij,  P^  and  | n^ |  are 

corrected  by  graphs  of  Figs.  15  and  16.  Pharos  of  and  >tj  weakly  depend  on 

conductivity  of  earth. 

Deviation  of  magnitudes  of  and  Pj  from  values  in  Figs.  17-19  due  t.c  change 
of  03  ./o)  and  t,  which  can  take  place  in  periods  of  ionospheric  perturbations 

*  3  * 

(flashes  on  sun),  is  determined  with  help  of  graphs  of  Figs.  12-14. 

Thus  with  help  of  basic  working  graphs  of  Figs.  17-22,  corrected  by  graphs  of 
Figs.  12-14,  allowing  for  deviations  of  air/u)  and  t,  graphs  of  Figs.  15-16  for 
calculation  of  conductivity  oQ,  formula  (4.4)  for  calculation  of  changes  of  earth 
radius  a,  and  formula  (2.98)  for  calculation  of  heights  of  points  of  radiation  and 
reception,  it  is  possible  to  calculate  all  characteristics  of  normal  waves  necessary 
for  calculation  of  distant  field. 


CHAPTER  V 


CALCULATION  OP  REGULAR  DEPENDENCIES  OP  AMPLITUDE  AND 
PHASE  OP  DISTANT  VLW  FIELD  ON  TIME  AND  DISTANCE 


§  1.  Calculation  Formulas  for  Components  of  Distant 
- VLW  Field  TMod'erTT - - - 


As  already  was  indicated  in  Chapter  I,  adequate  theory  of  [VLW]  (C.HB )  propagation 
can  be  formulated  only  on  the  basis  of  model  of  medium  considering  temporal  and 
horizontal  variations  of  parameters  of  medium.  Therefore  we  Bhall  use  subsequently 
calculation  formulas  for  model  C.  Using  formulas  (2.25) -(2.30) ,  with  function 
| B/A  |  determined  by  (2.128),  we  obtain  following  calculation  formulas  for  components 


of  fields 


-w-  »  _ 

(«*'»*) Dy(a’s) e  * 
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where  nj(0)  designated  coefficient  of  excitation  at  point  of  radiation  0,  and  nj(P) 
Is  same  at  point  of  reception  P. 

Integration  in  exponential  factors  is  carried  out  with  respect  to  the  shortest 
"beams,"  connecting  point  0  with  point  P. 

As  can  be  seen  from  Chapter  IV,  during  calculation  of  field  in  range  of 
lu-20  kilocycle,  it  is  necessary  to  take  into  account  only  three  normal  waves: 
quasi-TH^,  -THg,  and  -TE^.  Therefore  for  calculation  of  field  on  surface  of  earth 
formulas  (5.1)  take  form: 


where 


--'A  | 

r;r- 

sin* 

r 

e'4\THt 

i  ///, 

.  77:,|  (qv/m]. 

H 

T 
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(5.2) 
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:  /■//,  • 

(  l 

//T,|  (pv/m). 
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=  j'  ■ 
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Integration  in  (5.2)  and  (5.3)  is  carried  out  with  respect  to  the  shortest  arc 
oi'  great  circle  connecting  points  0  and  P,  since  distortion  of  "learns"  TK^,  TH^ ,  an  1 
TE„  is  small. 
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A. 


Replacement  of  true  trajectory  of  beam  by  geodesic  is  permissible  if  path  of 
beams  TH^  TH 2,  and  TE1  form  angle  of  more  them  10°  with  twilight  belt.  In  case  of 
paths  for  which  beams  slip  along  twilight  belt,  for  exact  calculations  it  is 
necessary  to  consider  distortion  of  beams  due  to  refraction  and  also  change  of 
coefficients  of  excitation  at  point  of  radiation  n(0)  and  at  point  of  reception 
n(P)  due  to  lens  effect. 

All  parameters  entering  in  formulas  (5.2) -(5.5)  can  be  determined  from  graphs 
given -^n  Chapter  IV  according  to  data  on  earth's  surface  and  lower  layers  of 

ionosphere  along  calculated  route.  Data  on  lower 
layers  of  ionosphere  are  determined  from  near  field  of 
very  long  waves  by  method  presented  in  Chapter  III. 
After  all  characteristics  of  normal  waves  on 


route  are  determined,  calculation  of  fields  Er,  E0, 

E  ,  Hr,  H0,  by  formulas  (5.2) -(5.3)  is  carried  out 


Fig.  23. 


/  by  analytic  or  graphic  method.  In  Fig.  23  is  depicted 

vector  diagram,  obtained  through  use  of  graphic  method. 

For  Er  basic  role  is  played  by  TH^  and  TH^  waves. 
With  increase  of  uistance  0,  all  vectors  of  waves 
revolve  counterclockwise  and  decrease  in  magnitude.  Thus  each  of  them  forms 
trajectories  in  the  form  of  twisted  spirals.  In  view  of  fact  that  >  a 2H, 
angular  velocity  of  rotation  of  TH^  waves  is  more  than  that  of  THg  waves.  Therefore, 
moving  in  revolving  system  of  coordinates,  together  with  vector  TH^,  we  shall  see 
rotation  of  vector  TH2  clockwise  with  angular  velocity  determined  by  difference 
alH  "  ^2H*  Here  total  amplitdue  and  its  phase  with  respect  to  wave  TH1  will 
oscillate. 

§  2.  Spatial  Dependencies  of  Amplitude  and  Phase 
of  Distant  Field  (Uniform  Routes) 

A.  VLW  Propagation  by  Lay 

We  shall  start  with  simplest  case  of  VLW  propagation,  when  route  passes  above 
marine  surface  and  1b  illuminated  evenly  by  sun.  In  this  case,  according  to  Chapter 
III,  following  parameters  of  medium  apply: 


a  — 6370  km,  3#—oc,  A,  >  70  km,  =  0,3-10** 

»  KM,  wJ.J 
—  8*10° 


,(•1  -  2-10’,  ht  s  75  km,  =  2- 10’*,  vi*>  =  10\ 

tlf  Ml 
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At  distance  of  more  than  2500  km,  for  calculation  with  accuracy  of  10$,  it  is 


possible  to  disregard  all  normal  waves  except  TH^.  Therefore  calculation  formula 
(5.1)  for  Ef  takes  simple  form: 


(w>->  r) 


(5.4) 


Field  E^,  owing  to  smallness  of  n1H  is  practically  absent.  If  radiating  Hertz 
doublet  and  receiving  point  are  on  surface  of  earth,  that  is,  b  =  a  and  r  =  a,  then, 
since  D(a',  a)  =  1,  formula  (5.4)  will  be  simplified  still  further: 


£,=  -iA, 


r~w  ,  » 
y  SP*’"' 


(’.«•  +  t). 


(5.5) 


where  9  is  angular  distance  between  points  0  and  P. 

In  expanded  form  instantaneous  value  of  E  will  be 


whence  amplitude  of  distant  field  |Er|  i3  equal  to 


(5.6) 


(5.7) 


In  Fig.  24  are  given  dependencies  of  |Er|  on  angular  distance  9  for  W  =  1 
kilowatt,  for  three  frequencies  f  -  12,  15*  and  19  kilocycles.  In  the  same  place 
are  given  for  comparison  results  of  measurements  of  daytime  field  from  round-the- 
world  expedition  of  1921-1922  under  summer  conditions  [43]. 

From  Chapter  III  it  follows  that  daytime  summer  conditions  differ  little  from 
daytime  conditions  in  period  of  equinox;  therefore  it  is  permissible  to  compare 
results  of  calculation  and  experiment.  Phase  of  distant  VLW  field  in  daytime,  as 
follows  from  (5.6),  grows  evenly  with  distance  for  a0  >  2000  km. 

On  distances  less  than  1500  km  it  is  necessary  to  consider  higher  types  of 
waves.  In  Fig.  25  is  presented  analytic  curve  of  |E  (0)|,  taking  into  account 
quasi-TH1,  quasi-THg,  and  higher  types  of  waves.  In  the  same  place  is  given 
experimental  dependence  of  |E  (0)|,  taken  with  help  of  aircraft  flying  on  route 
from  England  (340  km  from  Rugby)  to  Cairo  [45]  in  the  summer  of  1950.  Since 
results  of  work  [45]  are  presented  in  relative  magnitudes,  then  on  basis  of  good 
coincidence  of  theory  and  experiment  for  distant  field  (Fig.  24)  we  assumed  that 
both  curves  coincide  for  large  values  of  9.  Here  also  is  revealed  good  coincidence 
of  experiment  with  theory. 
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theoretical 


It  is  necessary  to  note  that  in  existing  literature  on  long  waves  are  widely 
applied  the  phenomenological  formulas  of  Osteen-Kogan,  Puller,  Ekspenshid,  and 
others.  All  of  them  have  following  structure: 

(5.8) 

where  D  =  a0  —  linear  distance, 

P  -  HI  —  moment  of  Hertz  doublet, 

1 

X  —  wavelength,  m  =  ^  -  1. 

First  factor  of  (5,8)  determines  field  of  dipole  in  free  space,  the  second 
considers  damping  in  medium.  If  this  formula  is  compared  to  (5.7)  and  we  consider 
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(5.9) 

(5.10) 


then  sharp  distinction  in  dependencies  of  first  and  secmd  factor  on  frequency  u> 
is  immediately  apparent.  According  to  Fig.  18,  depends  on  frequency,  and 
weakening  of  distant  field  with  frequency  (Fig.  24)  is  caused  by  fall  of  |nlh| 
(Fig.  19).  In  phenomenological  formulas  (5.8),  conversely,  pre-exponential  factor 


increases  with  frequency,  and  weakening  of  distant  field  with  frequency  should  be 
provide  by  second  factor.  Thus  slope  of  drop  of  phenomenological  curves  |Er(8)|  is 
increased  with  frequency.  From  Fig.  24  it  follows  that  slope  is  practically 
independent  of  frequency.  Therefore  phenomenological  formulas  give  grossly  incorrect 
quantitative  and  qualitative  results.  Unfortunately  they  occupy  central  place  in 
all  monographs  and  textbooks  on  propagation  of  radio  waves  in  sections  on  long  waves, 
misleading  broad  radiotechnical  circles. 

For  VLW  propagation  in  daytime  above  land  it  is  necessary  to  consider 
conductivity  of  earth  Og.  Formulas  for  calculation  of  field  above  land  surface  are 
not  changed,  and  in  earth  or  under  surface  of  water  will  be 

etc. ,  (5.11) 


where  Er(a,  6)  is  determined  by  formula  (5.6). 

If  route  is  nonuniform  with  respect  to  ionosphere  and  ground  conductivity,  then, 
according  to  (5.2)-(5.3),  it  is  necessary  to  replace  exponential  factor  by 

O 

integral  and  quantity  n^  by  product  n^O)  ‘n^P) ,  where  in  coefficients  n1H(0)  and 

n1H(P)  it  is  necessary  to  consider  ground  conductivity  near  radiator  and  receiver. 

Formula  for  E„  takes  form 
r 


«,*«>) 


Iff' 


(5.12) 


As  can  be  seen  "takeoff"  and  "landing"  sites  for  radio  waves  play  considerably 
smaller  role  here  the  route  itself.  In  this  respect  very  long  waves  differ 
essentially  from  medium  waves,  investigated  by  group  of  Soviet  scientists  under 
leadership  of  academicians  Mandelstam  and  Papaleksi.  It  is  necessary  to  note  also 
the  distinction  in  shore  refractions.  On  VLW  there  is  real  shore  refraction, 
nonvanishing  with  distance  from  shore,  but  it  is  extraordinarily  weak,  as  can  be 
seen  from  insignificant  increase  of  Aa^  during  transition  through  land-sea  boundary 
(Fig.  15). 

Influence  of  conductivity  of  land,  predicted  by  theory  within  limits  of  errors 
of  measurements,  agrees  with  experimental  investigations.  In  above  cited  work  [45] 
are  described  changes  of  daytime  field  strength  of  radio  station  in  San  Francisco 
(f  =  23  kilocycles),  taken  on  ship  proceeding  from  England  to  New  Zealand  through 
Danama  Canal.  Field  strength  was  increased  by  approximately  twice  during  transition 
from  radioroute  from  San  Francisco  to  Panama  Canal  (almost  entirely  land)  to 
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radioroute  from  San  Francisco  to  the  Pacific  Ocean  (15°  S,  100°  W)  (entirely  sea). 
From  Figs.  15  and  18  we  determine  by  means  of  extrapolation  that  for  sea  (o0  = 
=  cd)  is  equal  approximately  to  2.2.  For  land  of  average  conductivity,  taking  into 
account  frequency  correction,  -  0.8.  Hence  Er(Cyuja  )^Er(uope )  *  ?*  **  *s 

necessary  to  note  that  in  Fig.  8  of  work  [4J]  on  directions  of  routes  are  plotted 
field  strength  reduced  to  5000  km. 


B.  VLW  Propagation  at  Night 

Now  we  shall  consider  VLW  propagation  for  period  of  vernal  equinox  under 
conditions  when  entire  route  passes  through  ncnilluminated  region  above  marine 
surface.  Disregarding  irregular  variations  of  parameters  of  lower  layers  of 
ionosphere,  caused  by  turbulence  and  other  factors,  we  shall  consider  that  along 
entire  night  route  medium,  according  to  Chapter  III,  is  characterized  by  following 
parameters: 

g  —  6370  km,  ®®»5  KM, «J,|  —  O.S.IO11, 

*JJ|-8.10»,  Jk,  =  90  km,  «tt-l,5.10,,l  rfJJ-5.10*. 

•W-0.IO*. 

For  distances  of  more  than  2000  km,  for  calculation  of  field  strength  with 
accuracy  of  10$,  it  is  possible  to  be  limited  to  three  waves:  quasi-TH^  quasi-TH2, 
and  quasl-TE^. 

In  Figs.  26  and  27  are  given  amplitude  dependencies  on  distances  |Er(9)|  and 
dependence  of  phase  on  distance,  arg  E r(0),  for  case  b  =  a  and  r  ~  a  with  above- 
mentioned  parameters  of  medium.  Calculation  was  conducted  with  first  formula  of 
(5.2). 

Characteristics  of  normal  waves  and  nj  are  determined  on  graphs  of  Chapter 

IV. 

Letters  A,  B,  and  C  on  Figs.  26  and  27  denote  curved  pertaining  respectively  to 

frequencies  f  =  1J  kc,  16  kc,  and  20  kc. 

First  of  all,  let  us  consider  amplitude  dependencies  (Fig.  26).  Oscilliatory 
variations  of  amplitude  of  total  field  Er,  as  are  those  of  all  other  components, 

are  caused  by  interference  of  TH^  and  THr,  waves.  Quasi-TE^  wave  has  phase  velocity 

close  to  that  of  quasi-TH1  wave  and  considerably  smaller  amplitude.  Therefore  it 
gives  long-period  shallow  pulsations.  On  small  distances  TH2  wave  dominates;  on 
large  TH^  dominates.  Cause  of  this  is  large  coefficient  of  excitation  and  attenuation 
factor  of  TH2  wave. 
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Fig.  25. 

In  Fig.  28  are  given  distances  9min  and  6^  at  which  amplitude  of  Er  turns 
into  minimum  (solid  curves)  and  maximum  (dotted  curves).  Figures  on  the  right 
denote  reference  numbers  of  extrema.  In  interval  of  angular  distances  marked  by 
heavy  lines  TH.^  and  TH2  waves  are  approximately  equal;  therefore  low  points  of  |Er| 
attain  greatest  depth  here.  These  regions  9  we  call  Sm^n  mir 

As  can  be  seen  from  Fig.  28,  for  certain  frequencies  there  exist  two  deep 
adjacent  minima  of  approximately  identical  magnitude.  Dotted  curve  shows  dependence 
of  [ Er ]  on  frequency  at  0minrain*  From  Fig.  28  it  follows  also  that  with  increase 
of  frequency  <a,  region  0min  min  rapidly  departs  from  radiator  and  depth  of  low 
points  increases. 

In  upper  part  of  Fig.  28,  above  band  9  w  0 min  min'  (luas^"^i  wave  dominates;  in 
lower  part  quasi-THg  wave  dominates.  At  large  distances,  when  (6  -  0minmlrl)  is 
great,  amplitude  of  pulsations  of  curves  |Er|,  depicted  in  Fig.  26,  attenuates,  and 
l.ight  field  of  (Er|  may  be  calculated  by  monomial  formula  (5.6). 
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Mean  amplitude  of  |Er|  for  averaged  pulsations  at  small  distances  0  <  ^m-pn  tn 
is  determined  by  formula  (5.6),  where  as  basis  is  taken  TH^  wave.  Mean  amplitude 
of  |Er|  for  0  >  9min m1n  is  determined  by  formula  (5.6),  where  as  base  is  taken  TH^ 
wave.  Since  attenuation  factors  of  TH^  and  TH1  waves  are  different,  mean  amplitude 
| E r |  drops  for  small  values  of  0  faster  than  for  large  ones.  Therefore  monomial 
phenomenological  formulas  are  not  useful  for  description  of  space  regularities  of 
night  field  in  whole  interval  of  Q. 

Complicated  interference  character  of  night  field  very  much  hampered 
experimental  investigations  of  space  regularities  of  VLW.  Therefore  in  work  [45] 
are  given  only  quasi-minimum  and  quasi-maximum  values  of  |Er|  for  frequencies  of 
16-30  kilocycles,  plotted  in  Fig.  26  by  dotted  line. 

Despite  the  fact  that  mean  amplitude  of  night  field,  determined  by  TH2  wave  at 
distances  6  <  0min  and  by  TH^  wave  for  0  >  min*  ds  larger  than  day  values 

(Fig.  24)  in  regions  of  9  close  to  6mln  mIn>  real  night  field  can  be  considerably  less 
t.han  day.  Example  of  this  is  sharp  fall  of  distant  field  of  radio  station  GBR  (Rugby), 
f  =  16  kilocycles,  during  reception  in  Moscow.  From  Figs.  24  and  26  it  follows  that 

~  =  3,  which  coincides  well  with  experimental  swing  of  diurnal 

Lr(H0Ub  )  50 

variations  of  station  GBR  (§  3,  Chapter  V),  which  is  located  at  angular  distance  of 
0  =  °'4  =  Gminmin  fr°m  Moscow* 

Another  example  of  fall  of  night  field,  which  was  earlier  considered  an 
anomaly,  is  Kauka  (Hawaiian  Islands)  -  Tokyo  route,  with  extent  of  9  =  1  radian 
(frequency,  17.7  kilocycles).  From  Figs.  26  and  24  it  follows  that  when  0=1 


rn  linn 


,  jE^I  falls  in  region  0, 


„  where  S  40 

Jr(H0«b) 


minmin’  —  F7 - 7  ~  3T  =  which  coincides  with 


duirnal  swing  of  variations  of  |Er|  given  in  book  by  Acad.  A.  N.  Shchukin  [1], 

Thus  anomalous  daily  variation  of  distant  field  of  VLW  is  always  obtained  for 

regions  9  ~  emin  min  ’  w^ere  there  occurs  interchange  of  TH^  and  TH^  waves.  It  is 

interesting  to  note  that  in  range  of  10-20  kilocycle  ratio  E  (^eHL.)A:r(lioqi  )  in 

region  0  ~  ®mj_n  does  not  depend  on  frequency  and  equals  approximately  3. 

In  regions  0  <  and  G  >  ^min  min  day  field  is  less  than  night.  This 

case,  naturally,  is  observed  more  frequently  in  practice  of  long-distance 
communication  than  with  the  preceding.  Therefore  it  is  considered  normal.  Examples 


cf  normal  cases  are: 


1.  Route  from  Rocky  point  (United  States)  to  New  Southgate  (England),  distance 
0  =  0.85  radian  (frequency,  17.150  kilocycles).  On  it  is  observed  regular  increase 
of  night  field 


JO 


-0.5. 


From  Fig.  28  it  follows  that  at  f  =  17  kilocycle  in  region  9 
interference  maximum. 

From  Figs.  26  and  24  it  follows  that 


O.85  there  is 


2.  Route  from  Annapolis  (United  States)  to  Moscow,  distance  0  =  1.22  radian 
(frequency  f  =  19  kilocycles;  NSS) ;  on  it  in  periods  of  equinox  is  observed  regular 
increase  of  field  at  night  (§  5  of  this  Chapter) 


,  .  -  0,3240.35. 

From  Figs.  26  and  28  it  follows  that  9  =  1.22  for  f  =  19  kilocycles  is  between 
third  maximum  and  minimum  of  interference  curve  |Er|  and  ratio 


EjLm**) 

f<ONk) 


a  0,3. 


5.  Route  from  Naven  (Germany)  to  Tokyo,  distance  9  =  1.4  radian  (frequency 
f  =  16.5  kilocycles);  on  this  route  there  regularly  is  observed  5.5-6-fold  gain  of 
night  field  as  compared  to  day  [1], 

From  Fig.  28  one  may  see  that  for  f  =  16.5  kc,  9  =*  1.4  radian,  at  night  there 
is  a  fourth  interference  maximum.  According  to  Figs.  24  and  26 

j&SEli  a!?«I. 
fflMUllmf,  65  5* 

4.  Route  from  Hawaiian  Islands  to  Moscow,  distance  9  =  1.8  radian  (frequency 
17.5  kilocycles);  on  this  route  was  observed  4-fold  increase  of  night  field 
(authors'  measurement).  From  Fig.  28  it  follows  that  for  9  =  1.8  radian  and 
frequency  of  17.5  kilocycles  at  night  there  is  fifth  interference  maximum;  field 
in  it  is  of  the  order  of  50  microvolts/m  (Fig.  26).  By  day  (Fig.  24)  field 
constitutes  7  microvolts/m.  Thus 
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We  could  have  multiplied  number  of  examples  of  anomalous  and  normal  daily 
variation  of  |E  (t)|,  which,  within  limits  of  errors  of  experiment,  are  contained 
in  above-stated  theoretical  calculations. 

It  is  necessary  to  note  that  diurnal  swing  of  |Er|,  determined  by  ratio 

£ 

)  t  does  not  depend  on  conductivity  of  earth  along  route,  since  corrections 

r(H0Ub ) 

and  A6?H  are  identical  both  for  day  and  night.  Therefore  during  comparison 
of  experiment  with  theory,  we  used  data  for  oQ  *=  as. 

Now  let  us  turn  to  discussion  of  graph  of  night  phase  depicted  in  Pig.  27. 

In  it  are  given  phase  differences  between  total  vector  of  field  Er  and  vector  of 
basic  wave,  having  the  biggest  modulus  (Pig.  23).  Graphs  A,  B,  and  C  pertain 
respectively  to  frequencies  of  13,  16,  and  20  kilocycles.  Since  in  interval  of 
angular  distances  0.2  <  9  <  1.5  for  f  *  13  and  f  =  16  basic  wave,  according  to 
Figs.  26  and  28,  is  quasi-TH^  wa.e,  while  for  f  =  20  kilocycle  it  is  TH^  wave,  then 
for  graphs  A  and  B  on  ordinate  are  plotted  phase  differences  between  vector  of  total 
field  and  that  of  TH1  wave,  while  for  graph  C  ordinate  bears  vector  of  total  field 
and  of  TH2  wave. 

From  Fig.  27  one  may  see  that  deviations  of  phase  Acp  of  total  field  from  basic 
wave  oscillates  around  zero  value  and  attenuates  with  growth  of  8,  Maxima  of 
deviation  of  phase  Atp  from  zero  value  are  shifted  with  respect  to  maxima  of  amplitude 
'f  field  | E r |  by  1/4  period  of  space  variations,  that  is,  are  observed  at  various 
angular  distances.  With  increase  of  frequency  f  swing  of  oscillations  Atp  grows. 

At  our  disposal  there  are  no  experimental  data  on  measurement  of  phase  of 
distant  field.  Therefore  confidence  in  result.-  of  Fig.  27  can  settle  only  on 
coincidence  of  experimental  data  on  amplitude  of  field  with  calculated  values.  The 
latter  were  obtained  from  the  same  formulas  used  for  data  on  phase. 

§  3.  Temporal  Dependencies  (Diurnal  Variation)  of  Amplitude 
and  Phase  of  distant  Vlw  Field  on  Nonuniform  Routes 
of  Variable  Y1  lumlnanc"e 

For  practical  purposes  it  is  very  important  to  know  diurnal  variations  of 
amplitude  and  phase  of  distant  VLW  field.  In  this  section  we  shall  consider  two 
cases : 

1)  anomalous  dal  Ly  variation  :  rugby  (England)  -  Moscow  route  (r<  =  0.4,  f  = 

=  16  kilocycles), 


2)  normal  dally  variation  on  Annapolis  (United  States)  -  Moscow  route  (0  =  1.22 
f  =  19  kilocycles)  In  period  of  vernal  equinox. 

Regular  heterogeneities  along  route  are: 

1)  change  of  conductivity  of  earth  aQ, 

2)  changes  of  properties  of  ionosphere,  caused  by  shift  of  twilight  belt. 
Regular  change  of  phase  and  amplitude  of  distanl  VLW  field  at  fixed  point  of 

reception,  in  this  case  in  Moscow,  appears  only  because  )f  shift  of  twilight  belt, 
along  route. 

Calculation  of  field  of  Er  at  the  surface  of  earth  was  made  by  formula  (5.2). 
where,  instead  of  angular  distance,  is  introduced  average  time  on  route  t  by  the 
formula 


Here  0Q  is  angular  distance  between  0  and  P,  tQ  is  difference  in  local  times  between 
0  and  P.  Integration  is  carried  out  from  t  to  t  +  t&,  where  t  is  time  at  point  of 
reception  P. 

Dependencies  of  characteristics  of  normal  waves  Vj,  nj,  and  Hj  on  time  were 
determined  with  help  of  graph  of  heights  h^(t).  Fig.  9.  on  which  abscissa  t  was 
replaced  by  0  per  formula  (5.13)  and  from  Figs.  17-22  for  night  and  day  values  of 
Vj,  nj,  and  Kj.  Values  of  Vj,  *j,  and  in  twilight  period  were  obtained  from 
Figs.  17-22  by  means  of  linear  interpolation.  Since  curve  h^(t)  consists  of 
straight  lines,  then  integrals  entering  in  (5.2)  are  calculated  by  formulas  for 
areas  of  triangles. 

In  Fig.  29  are  given  computed  values  of  amplitude  and  phase  of  field  of  station 
GBR  (Rugby)  in  Moscow  for  period  of  vernal  equinox.  In  the  same  place  is  given 
experimental  curve  of  daily  variation  of  amplitude,  taken  on  19  March  1955  in  region 
of  Moscow.  All  data  on  amplitude  are  given  in  decibels  with  respect  to  daytime 
values  of  field,  which  is  distinguished  by  high  stability  of  amplitude. 

As  can  be  seen  from  Fig.  29.  there  is  good  coincidence  between  calculations 
and  experiment.  Theory,  in  virtue  of  simplified  Initial  data  (Fig.  9).  considers 
only  regular  oscillations  of  field.  Therefore  from  comparison  of  experimental  and 
calculated  curves  it  is  possible  to  dismember  regular  and  irregular  parts  of  field. 
This  will  be  done  in  following  section.  Here  we  are  interested  only  in  regular 
changes  of  daily  variation. 
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On  the  basis  of  presented  theory  it  is  simple  to  explain  temporal  variation  of 
amplitude  of  field  Er,  depicted  in  Fig.  29.  In  course  of  day  amplitude  Er  does  not 
change,  since  field  at  point  of  reception  is  caused  by  one  quasi-TH^  wave.  Small 
drop  of  amplitude  after  sunset  (point  A,  Fig.  29)  is  due  to  amplification  of 
adhesion  effect  for  wave,  due  to  rise  of  ionosphere  at  point  of  reception. 

This  drop  has  been  observed  by  a  number  of  authors,  but  until  now  has  not  been 
explained  [1,  2,  3,  4,  5]. 

With  advance  of  twilight  westward,  an  even  greater  part  of  route  i3  submerged 
in  darkness.  Therefore  field  of  Er»  determined  by  TH^  -wave,  is  strengthened,  since 
attenuation  of  TH.^  wave  at  night  is  less  than  by  day  (point  B) .  Only  after  twilight 
reaches  England  and  conditions  close  to  those  of  night  have  been  established  data 
point  of  radiation  does  THg  wave  start  to  reach  point  of  reception,  and,  owing  to 
mutual  extinguishing  of  TH^  and  TH2  waves,  total  field  drops.  The  biggest  drop  of 
Er  (point  C)  is  observed  under  conditions  for  which  TH^  and  THg  are  in  opposite 
phase.  In  period  of  ascent  of  sun,  along  route  the  picture  is  approximately 
reversed. 

Coincidence  of  daily  variation  for  calculated  and  experimental  amplitudes  forces 
us  to  assume  that  calculated  curve  of  dally  variation  of  phase  is  close  to  reality. 
The  biggest  deviation  of  phase  from  daytime  value  occurs  in  periods  of  sunsets  and 
sunrises,  when  vectors  of  and  th2  (Fig.  23)  are  perpendicular. 

In  Fig.  30  are  presented  calculated  curves  of  dally  variation  of  amplitude  and 
phase  of  field  of  station  NSS  (Annapolis)  as  received  in  Moscow  for  period  of 
vernal  equinox.  In  the  same  place  are  given  experimental  curves  for  amplitude 

t  : 

taken  in  Moscow  on  19  and  20  March  1955. 

In  daytime  at  point  of  reception  quasi-TH^  wave  dominates.  Evening  drop  at 
point  A  caused  by  adhesion  effect,  is  more  sharply  expressed  here  than  in  preceding 
case.  Amplification  of  field  at  sunrise  at  point  of  reception  and  after  sunset  at 
point  of  radiation  is  caused  by  interference  of  TH^  and  TBg  waves.  In  these  period 
there  occurs  exchange  of  leading  waves.  At  night  quasi-1H2  wave  dominates,  since 

0  <  0min  min  (Flg8*  26  28>* 

Divergence  of  calculated  and  empirical  curves,  as  follows  fronj  comparison  of 
two  empirical  curves  obtained  on  different  days  is  caused  by  irregular  processes 
in  ionosphere.  Diurnnl  swing  of  phase  of  NSS  is  more  than  4  times  more  than  in 
1  receding  case  (GBR),  which  is  caused  by  corresponding  increase  of  distance  and 
frequency. 
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§  4 .  Reception  of  VLW  Under  Water 


For  radiator  located  on  land  (b  =  a)  and  receiving  point  immersed  at  depth 
H  «  a  -  r  under  land  or  under  water ,  formulas  (5.1)  for  fields  taking  following 
simple  form: 


_  ,  L,  -»«i — 
£»!^'  *  * 

(5.14) 

• 

«,(«,*)  sfff, («,»)*  v  * 

(5.15) 

E(a» 

(5.16) 

(1  -/)// 

(5.17) 

Ef  [Hfi)  =>£f(a,  8).  e  • 

(5.18) 

Hq{H,*)=  HQ{afi)e 

(5.19) 

where  Er(a,  0),  H^(  a,  6),  Hr,(a,  9),  E^(a,  9),  and  H0(a,  0)  are  total  components  of 
electromagnetic  field  above  land  surface,  calculated  by  formulas  (5.2)  and  (5.3). 

X  -  wavelength  in  free  space  (in  meters), 

°0  -  conductivity  of  water  or  soil  ^^gteTT . 

With  error  of  the  order  of  factor  ,  it  is  possible  to  consider  that  in 

air 

-  £,-//,  •  120*;  Ef  -  Hf  •  120*. 

Therefore  on  graphs  (Figs.  26  and  ?7)  for  fields  of  Er(a,  0),  given  in  Chapter 
V,  it  is  easily  possible  to  determine  Er(H,  0),  H^H,  0),  and  Eq(H,  0)  in  earth  by 
changing  scales  on  ordinates  accordingly: 


for  E,, 

60Xsc 

for  H f  , 

(5.20) 

120« 

1  * 

'■"■r  E 6. 

1  —  * 

/6OXa0 


Remaining  components  in  this  work  are  no'  given,  and  according  to  the  formula 
(5.1),  depend  on  polarization  factor  jij . 

Since  for  sea  water  mean  value  Oq  ~  4  ^^y^eTeT*  c^en  for  mlddle  of  range  of 
VLW  (X  =  20,000m)  parameter  60 Xo^  =■  4.8«10^.  Therefore,  a-r'-’ording  to  (5.14)  and 

(5.16), 

E,(H,Q)  1 _ —  0  4fi  10“\  that  is 

component  E()  by  more  than  2000  times  exceeds  component  Er. 

A 

Since  does  not  drop  below  by  day  and  does  not  exceed  1/10  at  night, 

1  i 

component  constitues  E0  at  night  and  E 0  by  day,  that  is,  is  always 

greater  than  Er#  From  these  data  the  conclusion  can  be  made  that  the  most  effective 
underwater  and  underground  antennas  will  be  systems  reacting  to  component  E^, 
equivalent  to  component  H^. 

Emf  induced  by  distant  VLW  field  in  vertical  loop  directed  towards 
radiator  is  determined  by  formula 


where 


dHdx, 


10";  h/m. 


(5.21) 

(5.22) 


Integral  (5.21)  is  selected  with  respect  to  area  of  loop.  E^nd  Is  due  basically 
to  component  E^.  If  loop  has  rectangular  form,  with  side  l  parallel  to  surface  of 
earth  and  with  side  d  normal  to  it,  then 


( ! .23) 


where  Hq  is  vertical  coordinate  of  center  of  loop. 


With  loop  of  sufficiently  small  dimensions 


—IH.  —I  Hi 

En~-l»toHf{afi)Se  =£w(«,l)e  ,  (5.24) 

where  S  =  Id  is  area  of  loop, 

Eind(a,  0)  is  emf  in  loop  located  above  ground. 

Magnitude  of  6  for  X  =  20,000m  and  oQ  =  4  is  equal  ~0.5i  that  is,  damping  with 
submersion  of  loop  is  ~l/2  neper  per  meter  of  depth  or  4  db  per  meter.  In  loop 
located  in  vertical  plane  perpendicular  to  [beam]  direction,  in  radiator  will  be 
induced  emf  caused  by  component  (or  H0).  From  above-cited  data  it  follows  that 
Eind  of  such  loop  will  be  1/100  (by  day)  and  1/10  (at  night)  of  Elnd  of  loop  located 
parallel  to  beam. 

Significant  result  of  theory  expounded  here  is  fact  that  field  at  depth  H  under 
surface  of  water  (land)  depends  only  on  field  on  surface  of  water  (land)  at  point 
with  same  geographic  coordinates.  From  Chapter  V  it  follows  that  field  on  surface 
of  water  can  be  calculated  with  great  accuracy.  Therefore,  using  conversion 
formulas  (5.20),  we  determine  with  the  same  degree  of  accuracy  fields  under  water 
and  land. 

From  calculations  given  in  Chapter  V  it  follows  that  to  provide  reliable 
communication  on  VLW,  it  is  necessary  to  consider  regularly  sharp  drops  in  signal 
strength  observed  in  certain  intervals  of  distances,  depending  on  frequency  of 
signals  and  also  on  seasonal  and  liurna  L  variations  of  ionosphere. 
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Designations 


A 


a 


B 


B. 


B 

A 


—  scalar  electrical  potential, 

—  scalar  electrical  potential  of  spherical  k  layer, 

—  scalar  electrical  potential  <>“  normal  wave  of  number  ,j , 
=  0.1829- IQ'20; 

—  coefficient  of  focusing  of  beam, 

—  radius  of  earth, 

—  scalar  magnetic  potential, 

—  scalar  magnetic  potential  of  spherical  k  layer, 

—  scalar  magnetic  potential  of  normal  wave  of  number  J, 

—  potential  vector  function. 


—  potential  vector  function  for  spherical  k  layer. 


b  —  radial  coordinate  of  Hertz  doublet, 

10 

c  —  velocity  of  light  =  3.10  cm/sec, 

.  .  .  i  i . 

t  (x,  y)  ...  D  (x  ,  y  )  —  analytic  functions  of  two  arguments,  composed  of 
yj  j  cylindrical  Bessel  and  Neumann  functions  of  order  Vj 

and  their  derivatives. 


D0,  D1 ,  D2 
A,  D 

n 


r  BB  ' '  '  DTB 


lengths  of  paths  of  beam, 

distance;  name  of  lower  layer  of  ionosphere, 

diameter  of  heterogeneity  in  ionosphere, 

indices  of  refraction  and  conversion  of  beams  in 
ionosphere , 

matrices  of  refractive  index  in  1-t.h  and  2-nd  layer: 
of  ionosphere, 


n- 


d  —  distance  between  routes  of  waves, 

E  —  electric  field  strength, 

E  s  E  ;  Eg  —  components  of  electric  field  strength  in  spherical 
v  coordinates, 

Ir  U  Jr 

E„j  E  :  E„  —  components  of  electric  field  strength  in  layer  of 
r  9  9  number  k, 

Efnd  —  [emf]  (0.SC)  induced  in  loop, 

e  —  charge  of  electron, 

F  Q  —  operator  determining  field  according  to  data  on  medium 
r*  and  exciting  [factors]  of  number  k, 

f  —  frequency  of  oscillations, 

H  :  H.  ;  H„  —  component  of  magnetic  field  strength  in  spherical 
41  coordinates, 

Jr  Vr  Jr 

H:  H  :  H0  —  components  of  magnetic  field  strength  in  layer  of 
1  9  0  number  k, 

H0,  Hq  —  intensity  of  constant  magnetic  field  of  Earth, 

H  =  a  -  r  —  depth  of  submersion  of  point  of  reception  P  under 
surface  of  sea  water  or  land, 

Hq  —  depth  of  submersJ  on  of  center  of  receiving  loop  under 
surface  of  water  or  land, 

h  —  instantaneous  coordinate  of  height  of  ionosphere, 
measured  from  surface  of  earth, 

hl5  ^2  —  *ie^8h't  a^d  2-nd  layer  of  ionosphere,  measured 

from  surface  of  earth, 

K  —  averaged  height  of  turbulent  ionosphere, 

T  —  density  of  outside  current, 

Ir  —  radial  component  of  density  of  outside  current, 

I.  —  component  of  radial  component  of  density  of  outside 
J  current  I  in  normal  wave  of  number  J, 

Jm  —  imaginary  part  of  quantity, 

|5|  —  vector  function  of  current  density, 

Jv  —  modified  Bessel  function  of  order  v, 
j  —  reference  number  of  normal  wave, 
k  —  reference  number  of  spherical  layer  of  medium, 
kQ  —  wave  number  in  earth, 
k^  =  ~  —  wave  number  in  atmosphere, 

Lr  —  operator  of  normal  waves  on  coordinate  r, 

* 

Lr  —  conjugate  operator  of  normal  waves. 
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—  operator  or  wave  equation  on  coordinate  6, 

—  operator  of  wave  equation  on  coordinates  P  and  <p. 


L(i). 

s  • 


Legendre -Hansel  inunctions  of  order  of  first  and 
second  kind, 

—  dif  ferential-mr-'  rix  expression  generating  operator 


—  differential  expression  generating  operator  L^, 

—  constant  finite  quantity. 


—  mass  of  electron, 

—  concentration  of  electrons, 


Nj  —  normalizing  factor  for  normal  wave  >f  number  j, 

nj  —  coefficient  of  excitation  of  normal  wave  of  numier  j, 

n.„  —  coefficients  of  excitation  of  quasi-TE.  and  quasi-TH., 
waves,  J  3 

n  —  modified  Neumann  function  of  order  v, 

V 

0  —  point  of  radiation, 

P  —  point  of  reception, 

P  —  electrical  moment  of  Hertz  doublet, 

P  —  pressure, 

P  —  first  fundamental  function  of  Legendre's  equation 
vj  of  order  v,, 

tJ 

Q  —  second  fundamental  function  of  Legendre's  equation 
j  of  order  Vy 

Re  —  real  part  of  quantity, 

RpB  —  indices  of  reflection  and  conversion  during  reflection 
of  beams  from  layers  of  ionosphere, 

r  —  radial  coordinate, 

r„_  —  indices  of  reflection  of  beam  from  Earth, 
rn 

S  —  area  of  loop, 

Sj  —  function  of  elkonal  of  normal  wave, 

T  —  temperature, 


SJ" 


TE  —  waves  with  transverse  electrical  polarization, 

TH  —  wave  with  transverse  magnetic  polarization, 
t  —  time, 

U  —  component  of  eigenvector  of  conjugate  operator  of 
normal  waves, 

T>)  —  amplitude  of  modulated  normal  wave  of  nur.ter  j. 
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--  vector  function  of  conjugate  operator  of  normal  waves. 


V  - 


v  — 

vc  - 

W  - 

Yj  Z  - 


component  of  eigenvector  of  conjugate  operator  of 
normal  waves, 

velocity  of  electrons, 

velocity  of  propagation  of  signal, 

emissive  power. 


V  Zjk 


components  of  vector  function 

Y 


of  operator  of  normal 


a  — 

a  — 

a,  - 


“jH'  “ 


P  - 

- 


PJH;  PjE 

y 

A  and  5 

AJ 

n 


e  — 


e(k)  _(k) 
trr  '  e00  - 


e  - 


K.  - 


njH'  hje 


x  - 


^  - 


*  ,  i 

v  ,  V,  + 


yjH;  VJE 


components  of  eigenvector 
waves  in  k  layer, 
azimuthal  angle, 

angle  of  incidence  of  beam  on  Earth, 

real  part  of  wave  number  of  normal  wave  of  number  j, 

real  parts  of  wave  numbers  of  quasi-THj  and  quasi-TEj 
waves, 

angle  of  incidence  of  beam  on  ionosphere, 

attenuation  factor  of  normal  wave  of  number  j, 

attenuation  factors  of  quasi-THj  and  quasi-TEj  waves, 

angle  of  slip  of  beam  upon  striking  ionosphere, 

increments  and  errors, 

Wronskian  for  Legendre's  equation, 

attenuation  factor  of  radio  waves  in  land  anti  in 
water, 

complex  dielectric  constant, 

complex  dielectric  constant  of  Earth, 

components  of  tensor  of  complex  dielectric  constant 
in  layer  of  number  k, 

angular  spherical  coordinate, 

polarization  factor  of  normal  wave  of  number  j, 

polarization  factors  of  quasi-THj  and  quasi-TEj 
normal  waves, 

wavelength, 

eigenvalue  of  conjugate  operator, 
small  parameter, 

complex  wave  number  of  normal  wave  of  number  j,. 

complex  wave  numbers  of  quasi-TH.  and  quasi-TE. 

J  J 

normal  waves. 
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veff  —  number  of>  collisions  for  electrons, 

veff  —  nuraber  °?  r*ollisions  in  k  layer, 
f-  =  cos  9, 

t r  =  J.1416, 

£  —  summation  operator*-* 
a  —  specific  elect  rical  conductivity, 

Oq  —  conductivity  of  earth, 
i  —  time  or  propagation  time  of  signal, 
t  —  parameter  of  anisotropy  of  ionosphert , 

4>j  —  phase  factor  of  forced  normal  waves, 

<p  —  angular  spherical  coordinate, 

<p  —  phase  angle, 

Xj  —  eigei. value  of  number  J  operator  of  normal  waves, 

—  phase  factor  of  free  normal  waves, 
a)  —  angular  frequency  of  oscillations, 

oo0  —  critical  frequency, 

<a0k  —  critical  frequency  of  k  layer, 

—  gyromagnetic  frequency, 

“r 

— - parameter  of  reflectance  of  ionosphere. 

A.  y.  —  division  of  goniometer  (0.001  of  distance). 
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